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ABSTRACT 



For (planar) closed self-avoiding loops we construct a "holographic" map 
from the loop equations of large- N QCD 4 to an effective action defined over 
infinite rank Hitchin bundles. We compute the effective action in a closed 
form in terms of functional determinants. The construction works as follows. 
Large- iV QCD 4 on T 2 x T 2 (T 2 x Tq) is equivalent by Eguchi-Kawai to a par- 
tially quenched (twisted) theory on T 2 described non-canonically by Hitchin 
pairs (A = A z dz + Agdz, \I/ = —iD u dz + iD u dz) with A transforming as a 
connection and $ as a Higgs field with respect to the gauge group on T 2 . 
(Centrally extended) parabolic Hitchin bundles classify a dense set of orbits 
of pairs (A, \1/ = ift + ■0) obtained as hyper-Kahler quotient by the action 
of the gauge group. We study the loop equations of the quenched (twisted) 
theory for the non-Hermitian connection B = A + i^Ktp + kt/>); k / and 
planar loops. We change variables in the loop equations from (centrally ex- 
tended) parabolic Hitchin bundles to the corresponding (centrally extended) 
holomorphic (i.e. chiral) de Rham local systems. The key point is that the 



quantum contribution to the loop equations for the holomorphic de Rham lo- 
cal systems can be evaluated by residues in a certain regularization and turns 
out to depend on the loop orientation and self-intersection. The outcome is 
that, for self- avoiding loops, the original loop equations are implied by the 
critical equation of an effective action (the master equation) defined over the 
microcanonical ensemble of the inductive limit of finite dimensional gauge or- 
bits in the quenched theory, or over infinite dimensional orbits, arising from a 
central extension of Hitchin bundles, in the twisted one. Finally, on the basis 
of the geometric Langlands correspondence, we propose a dual construction 
in which the moduli of Hitchin bundles are quantised and become differen- 
tial operators. The loop equations are then implied by the condition that 
the master equation annihilates the vacuum and the master field is realized 
as an element of a quasi-conformal (vertex) algebra. 



1 Introduction. 



To say it in a nutshell, in this paper we construct an effective action whose 
critical points solve the loop equations [f , 2] of large- N QCD in four dimen- 
sions [3] for (planar) closed self-avoiding loops. 

Being a functional over Hitchin systems [4, 5] (see appendix B), our effective 
action is, in a sense, the analogue for large- N QCD 4 of the Seiberg-Witten 
low energy effective action for the Abelian branch of Af = 2 SUSY gauge 
theories in four dimensions [6]. However, we should stress that our effective 
action is derived directly by the loop equations of the large- N pure gauge 
theory in four dimensions for (planar) closed self-avoiding loops, once care- 
fully chosen variables are employed in the functional integral. 
The effective action is derived via a "holographic" map [7, 8, 9], that allows 
us to map the loop equations, that live on the boundary, i.e. on the loops, 
into a master equation, that is the critical equation of an effective action 
defined on the bulk, i.e. on the space-time. 

One key feature responsible for the existence of the "holographic" map is a 
" crural" change of variables, that allows us to reduce the computation of the 
quantum contribution in the loop equations to the evaluation of a residue, 
that depends only on the loop orientation for self-avoiding loops. 
Another important feature of our "holographic" map is that the residue, that 
arises from the contour integral along planar self-avoiding loops in the loop 
equations, is cancelled by a u(l) anomalous contribution from the "chiral" 
change of variables in the bulk, provided the conformal factor of the metric 
is suitably chosen. 

Despite the holographic suggestions, in this introduction we would like to 
make contact with the concept of master field [10], that would occur as the 
solution of the master equation for our effective action, with special regard 
to some algebraic features that arose in the subject. 

From a purely algebraic point of view, it is of the utmost importance that the 
algebra generated by the monodromy of the connection that solves the loop 
equations for closed self-avoiding loops be hyper-finite, i.e. be the weak limit 
of a sequence of finite dimensional von Neumann algebras (see appendix A), 
a property shared (trivially) by the algebra generated by the distribution of 
the eigenvalues that solves the large- N limit of some matrix models [11]. For 
the existence of an effective action defined over a hyper-finite algebra (a HF 
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effective action in short) a crucial role is played by the condition that closed 
loops be self-avoiding in the loop equations. Indeed it has been known for 
some time that in the general case of arbitrary self-intersection, the ambi- 
ent algebra, in which the connection that solves the loop equations lives, is 
(algebraically isomorphic to) a free group factor with an infinite number of 
generators (see appendix A), that is not a hyper-finite von Neumann alge- 
bra. The lack of hyper-finiteness seems to be one of the reasons that so far 
prevented an exact treatment of the general case. 

The existence in large- N QCD4 of an hyper- finite effective action for (pla- 
nar) closed self-avoiding loops, as opposed to the general case, has a precise 
analogue in large- N QCD2 that is explicitly solvable [12]. In the two di- 
mensional case the von Neumann algebra generated by non-overlapping but 
self-intersecting loops at a point is precisely the free group factor with an 
infinite number of generators [13, 14, 15, 16], a type 1I\ non-hyper finite von 
Neumann algebra [16]. However functional techniques show that the master 
field, that solves the loop equations of large- N QCD2 for closed self- avoiding 
loops, lives instead inside a commutative hyper-finite von Neumann algebra 
and indeed it can be determined by the distribution of the eigenvalues [17]. 
This fact suggested to the author to look for hyper-finite solutions of the 
loop equations for closed self-avoiding loops in the four dimensional case. In 
a sense this approach follows the original Witten proposal for the master field 
[10], that contained implicitly some kind of finite dimensional approximation 
property in the prescription of evaluating gauge invariants operators on the 
master field by means of the normalised approximately finite dimensional 
matrix trace liniTv -^Tr^. In fact the existence of an approximately finite 
dimensional trace is weaker than hyper-finiteness but implies that the von 
Neumann algebra of the corresponding GNS representation is embeddable in 
the factor [18], the ultra power of the unique hyper-finite factor R (see 
appendix A). 

The plan of the paper is the following one. Sect. (2) is an algebraic detour in 
which we recall some old and new ideas about the master field. Though it 
is not strictly necessary for the rest of the paper it helps to put sect. (3) and 
sect. (4) in the proper algebraic perspective. We make also some comment on 
the physical meaning of restricting to closed self-avoiding loops. In addition 
in sect. (2) we recall the quenched [19, 20] and twisted [21] Eguchi-Kawai 
(EK) [19] reductions that are known to be equivalent to the large- N limit of 
QCD. The quenched and twisted theories that we refer to are the continuum 
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version of their original lattice formulation [19, 22, 23]. A partially quenched 
[24] or twisted theory [21] is used in this paper as a technical tool to embed 
an inductive sequence of finite dimensional Hitchin systems into the QCD^ 
functional integral. 

Sect. (3) contains a first change of variables, from the gauge connection to 
the moduli fields of finite dimensional gauge orbits obtained by hyper-Kahler 
[4, 25] reduction in a partially quenched Eguchi-Kawai model or to the moduli 
fields of the hyper-Kahler reduction of a certain infinite dimensional central 
extension of Hitchin systems in the twisted case. In this section the func- 
tional integral is represented as an integral over the (energy) levels of the 
microcanonical ensemble for the action of the local gauge group. However 
an effective action on the gauge orbits, that turn out to be Hitchin systems 
(or an infinite dimensional central extension of, in the twisted case), cannot 
be derived by the loop equations for these variables yet. 
In sect. (4) the effective action is constructed by means of a further change 
of variables, from the moduli fields of the orbits of the SU(N) gauge group 
to the moduli fields of the corresponding holomorphic (see appendix B) de 
Rham local system (or a central extension of). 

Objects related to the holomorphic de Rham local system have been intro- 
duced under the name of opers by Beilinson and Drinfeld in their study of the 
geometric Langlands correspondence [26] (opers will play a role in sect. (6), af- 
ter quantising Hitchin systems as an alternative tool to solve the loop equa- 
tions). The construction of the effective action involves a de Rham local 
system for a partially quenched or twisted Eguchi-Kawai model. Hyper- 
finiteness is shown essentially by construction in the quenched version, while 
it follows with little more effort in the twisted theory, since in this case 
the gauge algebra has to be defined from start on an infinite dimensional 
Hilbert space and only later on approximated by compressions that define 
finite dimensional representations of the fundamental group. In any case 
hyper-finiteness should hold also in the twisted theory, at least because of 
the equivalence of the quenched and twisted theories in the large- N limit. 
In sect. (5) the effective action is computed in terms of functional determi- 
nants. The effective action is the sum of three terms. The first one is the 
classical action of the EK partially quenched or twisted model. The second 
one is the logarithm of a functional determinant that arises from localisa- 
tion on the microcanonical ensemble. The third one is the logarithm of an 
" anomalous" Jacobian that arises from the change of variables from the mod- 
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uli of Hitchin systems to the moduli of the holomorphic (i.e. chiral) de Rham 
local systems. Despite the techniques employed in this paper are mostly 
based on two dimensional differential geometry, we find that the localisation 
determinant on the microcanonical ensemble of the Hitchin systems can be 
computed in a manifestly four dimensional Euclidean (or Lorentz) invariant 
way. Remarkably we find that the first coefficient of the beta function is 
exactly reproduced, at lowest order in the expansion in powers of the moduli 
fields (that at lowest order coincides with the usual perturbative expansion), 
by the contribution of the localisation determinant in the master equation. 
It follows immediately that the techniques employed here apply also to the 
M = 1 super-symmetric extension of large- N QCD4 obtained adding to the 
pure gauge theory minimally coupled Majorana gluinos in the adjoint rep- 
resentation, since in this case the N = 1 effective action is obtained simply 
adding the logarithm of the fermion Pfaffian to the effective action of the 
pure gauge theory. 

In sect. (6), using ideas based on the geometric Langlands correspondence 
[26, 27, 28, 29, 30, 31, 32] (and references therein), we suggest another way 
of solving loop equations for closed self-avoiding loops by quantising Hitchin 
systems. In this last formulation, due to quantisation, hyper-finiteness of the 
ambient algebra may get lost, keeping nevertheless (quantum) integrability. 
In this case the master field becomes an element of a quasi-conformal (ver- 
tex) algebra. This is very close to interpreting the master field as an element 
of a vertex algebra of a conformal field theory, as it would be required by the 
conjectured gauge- field/string duality [33, 34], but for the fact that our alge- 
bra is obtained first solving the Hitchin self-duality equations on a punctured 
torus and then quantising the resulting moduli space, while in the usual con- 
struction of conformal blocks first the space of two dimensional connections 
is quantised and then the quantised Hitchin equations are imposed as con- 
straints on the space of physical states [35, 36]. 
Our conclusions are contained in sect. (7). 

In appendix A we recall some elementary facts about von Neumann alge- 
bras, that may help to understand the revival of the algebraic approach to 
the master field. 

In appendix B we remind the reader of several notions of geometrical nature, 
about Hitchin bundles, that are used in the paper. 

In appendix C we compare our approach to the loop equations with Witten 
localisation formulas in QCD 2 [37, 38, 39, 40] . 
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In appendix D we compute, in a manifestly Euclidean (or Lorentz) invariant 
way, the divergences of the localisation determinant at lowest order in the 
expansion in powers of the moduli fields. 



2 Loop equations: an algebraic detour. 



Recently there has been a revival of applications of ideas in the theory of von 
Neumann algebras (for a comprehensive review see [41]) to field theory, in 
particular to the large- N limit (see for a recent review [42]). We regard the 
large- N limit as defined by the loop equations [1, 2], that are the Schwinger- 
Dyson equations in the large-iV limit for the expectation values of traces 
of Wilson loops in gauge theories or for a complete set of gauge invariant 
operators in other theories (i.e. matrix models). They are obtained noticing 
that the integral of a derivative vanishes: 



where the sum over the index a is understood. Here A a are Hermitian gen- 
erators of the Lie algebra and 



ty(x,x;A) is the monodromy matrix of the connection along the closed 
loop C(x, x) based at the point x, i.e. ^(x, x; A) is the Wilson loop. For the 
group SU(N) using the identity: 




(1) 






(3) 



we get: 



= /D^exp(--/ 



TrF^ u d 4 x)(—Tr(D fM F fMU (z)^(x, x; A)) + 



5 



+i 



J C(x,x) 



f dy u 5 {d \z - y)Tr(V(x, y; A))Tr(V(y, x; A)) + 



-±-Tr(*(x,y;A)*(y,x;A))) (4) 



where the last term may be disregarded in the large- A" limit. The second 
term is the quantum contribution to loop equations, while the first term is 
the classical one. Thus, using the factorisation of gauge invariant operators 
in the large- iV limit and noticing that the expectation value can be combined 
with the matrix trace to define a new generalised trace r, our problem is to 
find a (unique) solution A^x) to 



for every closed contour C, with values in a certain operator algebra with nor- 
malised (t(1) = 1) trace r [13, 14, 15]. Such a solution is named the master 
field [10]. The recent revival concerns the precise identification of the alge- 
bra and the trace relevant for the large- N limit of field theories [13, 14, 15]. 
Operator algebras have already played some role in the early days of construc- 
tive field theory, but they have been later abandoned in favour of functional 
integral techniques as soon as it became clear that the algebra of a single 
Euclidean scalar field is commutative. In our case however the algebra of 
Euclidean matrix fields in the large- iV limit is highly non commutative. 
The recent interest in the interplay between the theory of von Neumann alge- 
bras and the large- iV limit was driven on the mathematical side by Voiculescu 
use [16] of large- iV matrix models as a tool in the study of non hyper-finite 
von Neumann algebras, in particular of the "elusive" free group factors that 
are left out by Connes [41] algebraic classification of all hyper-finite factors 
(see appendix A). On the physical side, after Singer suggestion [44] that type 
Hi representations of Voiculescu free group factors are the relevant ones for 
the study of the large- A" limit, several old and new results about the master 
field have been interpreted in the light of the theory of von Neumann algebras 
[13, 14, 15]. These investigations culminated with the precise identification 
of the ambient algebra and the trace for all matrix models in the large-A" 
limit, at least to all orders in the planar perturbation theory [13, 14, 15]. 
In this section we recall some old and recent ideas about the master field, 




= -t(D ii F iu ,(z)*(x,x;A)) + 
dyj {d \z - v)t(*(x, y; A))r^f(y, x; A)) 



(5) 
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interpreting them from the more recent algebraic point of view (see [42] for 
a review). Though this is possibly historically incorrect, it is useful for our 
purposes. After the discovery of large- N factorisation [43] of gauge invariant 
local operators, Witten suggested that the large- iV limit could be understood 
as localisation of the functional integral on a unique gauge orbit, the master 
orbit, on which gauge invariant operators should be evaluated by means of 
the usual normalised matrix trace in the large- iV limit [10] . In the retrospect 
we observe that the use of this matrix trace is tantamount to postulating that 
the von Neumann algebra generated by the master field is of type ll\ [44] and 
admits an approximately finite dimensional trace, in particular is embeddable 
into R", the ultra product of the unique hyper-finite factor R. If in addition 
the trace is uniformly approximately finite dimensional the associated GNS 
representation is hyper-finite (HF) (see appendix A) [18]. While the aim of 
this paper is to find observables in large- N QCD^ for which the original Wit- 
ten idea of a possibly hyper-finite master field can be concretely realized, in 
general the hyper-finiteness condition turned out too much restrictive, even 
in the case of the large- N limit of a free scalar field [45] . Haan [46] discovered 
that the von Neumann algebra of the master field for a free scalar field was 
the Cuntz algebra with an infinite number of generators (though he did not 
give it this name). Haan found the following representation for the master 
field of a free scalar field in the large- N limit: 



where | Q) is the vacuum state. In addition he found, on the basis of the 
large-iV Schwinger-Dyson equations, that also in the interacting case of a 
scalar theory with quartic interaction (but in a way that trivially extends 
to any interaction) the ambient algebra of the master field was the vacuum 
(i.e. tracial) GNS representation of the Cuntz algebra. The Cuntz algebra 
in general is generated by the following creation and annihilation operators, 
that satisfy: 



(f>(x) = J a(p)exp(ipx) + a(p)* exp(— ipx)dp 

with the operators a(p) and a(p)* satisfying: 

a(p)a(qy = 5 (d \p-q) 
a(p) | ft) =0 



(6) 



(7) 



did* = dijl 



(8) 
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and no other relation but the completeness condition: 



^2a*a { = 1 - Pq 



(9) 



where Pn is the orthogonal projector into the vacuum. A tracial representa- 
tion is obtained requiring: 



where | Q) is the vacuum. After Voiculescu work we know that the alge- 
bra generated by the self-adjoint elements a* + a { is a type ll\ algebra not 
hyper-finite for % > 2. The Cuntz algebra with more than one self-adjoint 
generator in its tracial representation is algebraically isomorphic to the free 
group factor with an equal number of generators, that was known already 
to von Neumann to be an example of a non-HF factor [41]. The occurrence 
of free group factors in the large-iV limit of field theories is implicit in the 
work of [47] on the effective action of the large- iV limit via the introduction of 
symbols that " neither commute nor anti-commute" (an informal definition of 
a free algebra) and in a way or another enters (implicitly) most of the subse- 
quent developments (for example in the use of the loop algebra [48, 49]) until 
it has been made explicit in the work of [13, 14, 15] through the introduction 
of the Cuntz algebra as the ambient algebra for the master field. 
In particular in [13, 14, 15] it is proved, to every order in planar perturba- 
tion theory, that the master algebra is a sub-algebra of the Cuntz algebra 
generated by the following elements: 



where M... are the planar correlation functions, so that the knowledge of the 
master field is in fact equivalent to the knowledge of the solution of the planar 
limit. To summarise, we know that the master algebra of a large- N free field 
is a non-HF Cuntz algebra and that the master algebra of an interacting 
large- N matrix model is a sub-algebra of the Cuntz algebra, probably non- 
HF too. However in some very special cases, such as the one-matrix model 
[11], traces of local operators can be reconstructed by of the distribution of 
the eigenvalues, that generates an Abelian HF algebra (in fact in this case 
the master algebra is the commutative von Neumann algebra generated by 
M M for just one y). In this case HF does not come as a surprise, since a 



di | Q) = 



(10) 



= a* -\- M^ v a u + M^ p a v a p + ... 



(11) 
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commutative C*-algebra is of type I and therefore HF (see Appendix A). 
We regard the lack of HF in the general case as the most disappointing, 
making difficult any progress in finding exact solutions of the loop equations 
(unless, perhaps, some kind of exact integrability is present). In this paper 
we renounce to attempt a solution of the large- N limit for all observables, but 
rather we concentrate on those physically relevant observables for which a HF 
solution may exist. In this search we use as a guide the important example 
of large- N QCD 2 that is explicitly solvable [12]. For QCD 2 a gauge exists in 
which the theory is Gaussian, so that the master algebra is generated by the 
Cuntz algebra in that gauge [13, 14, 15]. In addition the algebra generated 
by Wilson loops based at a point is the free group factor with an infinite 
number of generators, one for each simple loop, i.e. non overlapping self- 
intersecting loop at a base point [14]. This can be proved noticing [14] that 
the iterative solution of the Migdal-Makeenko equations with respect to the 
number of self-intersections coincides with the iterative equation associated 
to a Voiculescu free family of non- commutative random variables [16], i.e.: 

{(W(Ai) ~ (W(A 1 )))...(W(A n ) - (W(A n )))) = (12) 

This provides iterative equations for the expectation values of products of 
n random variables in terms of expectation values of a lower number of 
variables. The entire family is thus determined by the expectation values of 
self-avoiding loops. For each simple loop of area A& it is possible to give an 
explicit unitary master field in terms of the Cuntz algebra [14]: 

W(A k ) = exp(iy/A k (a k + a*)) (13) 

However the expectation value of a closed self-avoiding Wilson loop on a 
large sphere can be computed by means of the distribution of the eigenvalues 
and as such admits a master field that generates a commutative HF algebra, 
as the master field of the one-matrix model. Indeed, on a sphere of large 
area, the average of a Wilson loop of area A is given by [50]: 

r N 

w(A) = z~ 1 JU do, n(ft - e 3 ) ex P (-— £ ej) x 

x N^cosOi (14) 
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where: 

A 

Vi - Oj) exp(- 



r N 

z = jl[de l H(e % - ex P (-— £ $1) (is) 



In the large- N limit the distribution of the eigenvalues is determined by the 
saddle-point equation: 

In sect. (4) we will see that the existence of HF solutions of the loop equations 
is somehow linked to the condition of non self-intersection. But even in this 
case, considerable difficulties arise when a "simple" (i.e. not living in a free 
group algebra) solution of the loop equations is attempted. Indeed for a 
closed self-avoiding loop based at x the loop equations degenerate into: 

r N r N 

= J 2M„exp(-— J TrF^d 4 x)(-Tr((D,F, u (x) + 

+i f dy v 6W(x - y)Tr(l))*(x,x;A))) (17) 

JC(x,x) 

We see therefore that, despite the equations have been "linearised" in such a 
way that the classical and quantum terms have the same operator structure, 
the quantum contribution has still a different functional dependence on the 
loop C from the classical term, because the distribution: 

if dyj {d \x-y) (18) 

JC(x,x) 

is obviously loop dependent. This is even more manifest writing the loop 
equations for self-avoiding loops in the Migdal-Makeenko form [1, 2]: 

= g(x) S(Trn r C )) + 
da^(x) 

+ f dyJ {d \x-y)Tr(l)(Tr^(x,x;C)) (19) 

JC(x,x) 

Somehow anticipating the result of sect. (4) we will manage to change vari- 
ables in the loop equations in such a way that in the new variables the 
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corresponding distribution in the quantum contribution will be loop inde- 
pendent for self-avoiding loops. This will open the way to find an equivalent 
effective action. 

Since we want to solve loop equations only for closed self-avoiding loops we 
wonder how much physical information is lost by this restriction. Obviously 
we can reconstruct by such loops the string tension at large distances and the 
running of the coupling constant at short distances. The next problem is the 
spectral information, i.e. the meson or the glueball spectrum. This turns out 
to be dependent on the dimension of space-time. Essentially, for the meson 
spectrum, we want to compute the Dirac propagator in the background of the 
master field. By employing a path integral representation for the propagator, 
we represent it as an average over Wilson loops computed along Brownian 
paths (for a recent use of this representation in the framework of the large- TV 
limit of QCD see [51, 52]). Now in d = 2 Brownian paths self-intersect any 
number of times with probability one, while in d — 4 they are self-avoiding 
with probability one (see for example [53]). Therefore in d — 4 the spectral 
information is contained in self-avoiding loops. An analogous argument holds 
for the glueball spectrum. We should mention that the construction of the 
effective action that we present in this paper holds only for planar loops, 
while in general to reconstruct the meson or the glueball spectrum we would 
need also non-planar ones. Since, however, every loop can be made to lie 
on a surface, we think that our construction can be extended to the general 
case at the price of making the geometry more involved, though this will not 
be attempted in this paper. We end this section recalling the EK reduction 
[19], that plays a significant role in our approach to the loop equations. The 
quenched EK reduction [19, 20] is essentially a way of constructing a theory 
with the same master field, i.e. the same large- iV loop equations, but with 
fluctuations reduced by a factor of N" 1 . In the process space-time degrees 
of freedom are reduced to lower, eventually zero dimensions. The partition 
function of the quenched EK model is given by: 

Z = Jim f dA, exp(^^Tr [p, + A?,p v + Aff) (20) 

where the integration is now over constant N x N matrices and p^ are di- 
agonal matrices with uniform distribution of eigenvalues, i.e. p^ are the 
quenched momenta. In the Eguchi-Kawai quenched prescription the volume 
of the phase space divided by (27r) d is assumed to scale as iV according to 
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(2n)~ d A d V = N in such a way that the action is effectively rescaled by a 
factor of N in order to compensate for a reduction of a factor of N in the 
entropy of the functional integration, that now is only over constant matri- 
ces. One way of seeing that the EK model is equivalent to the large- N limit 
of QCD is to look at its loop equations: 

f N(2n) d 
= J dA^exp( 4 V AA j Tr\ Vil + A^p v + A u f) x 

x (jTr[ Pfl + Af, [p, + Af,p v + A v ]]*( x , x; A) + 
r A d 

_< L dy„-—jTr(*(x,y-,A))Tr(*(y,x;A))) (21) 
JC(x,x) {2ir) a 

Provided that the coupling constant was rescaled as in Eq.(20), these are 
the same loop equations as for large-^ continuum QCD, but for the fact 
that there is no delta function in the quantum contribution, since in the 
reduced model we are integrating over constant matrices. Therefore the loop 
equations of the reduced EK model coincide with the ones of the un-reduced 
theory provided the trace of open loop vanishes, that is a way of reproducing 
the delta function. This is ensured by the global R d symmetry (if it remains 
unbroken) (see [42]): 

Ap^ppl + Ap (22) 

that implies: 

jj(Tr(*(x, y; A))) = exp(ip(x - y))(^Tr(*(x, y; A))) (23) 

In fact the quenching prescription was introduced precisely for this aim. A 
variant of the quenched EK model is the twisted model [21]: 

r N 

Z = J dA^ exp(— Tr[T, + A^ T v + A v f) (24) 

where V ^ are (unbounded) self-adjoint operators acting on an infinite dimen- 
sional Hilbert space satisfying: 

[EV,r„] =i(r%,i (25) 
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with (6 1 ) fiu vanishing in the large- N limit. In modern language the twisted 
model is a gauge theory defined over the non- commutative Euclidean space: 

[x^x v ] = iO^l (26) 

The non-commutative Yang-Mills theory (for a review see [54]) equivalent to 
the twisted model is obtained working with commuting coordinates and re- 
placing the ordinary product between functions with the Moyal star product 
defined by: 

d d 

f(x) -k g(x) = exp(i9 lxl/ ——)f(x + y)g(x + z) \ y=z =o (27) 

Thus the partition function of the twisted theory can be re-casted into its 
non-commutative form, that looks like a matrix model: 

r N 

Z = JdA, ex P (-—Tr(t[C„ C u ] - (6~ l ), v f) (28) 

where C M = {—i0^ 1 )^ lu x L/ + iA^ and the product is the Moyal product. The 
loop equations of the non- commutative theory [55, 56] differ from those of 
the commutative theory for the need of taking an integral over the base-point 
of the Wilson loop, in order to get gauge invariant quantities, since the trace 
now includes integration over the base point of the Wilson loop. However, 
after taking the expectation value, the trace of Wilson loops is proportional 
to the volume of space-time, because of translational invariance, so that we 
get a normalised trace dividing by the volume. In the 9^ — > oo limit the loop 
equations of the non-commutative theory (for any positive N, even N — 1) 
coincide with the ones of the large- iV commutative theory. 



3 The microcanonical ensemble. 

It is quite well known that every Euclidean field theory can be thought as 
the partition function of a system of classical statistical mechanics in which 
the coupling constant squared is identified with the classical temperature 
(see for example [57]). Following this interpretation the Euclidean functional 
integral: 

/ exp(-\s(A))DA (29) 
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looks like the classical partition function: 

Jex P (-f3H(p,q))dpdq. (30) 

The second integral has an obvious representation as an integral over the 
energy levels of the microcanonical ensemble (see for example [58]): 

J exp(—PH(p, q))dpdq = 
= J exp(-f3E)5(H(p,q) - E)dpdqdE = 

= J exp(-(3E)Det(VHy 1 dE (31) 

The aim of this section is to find an analogous representation for the partition 
function of QCD^. The motivation is the following one. The original Witten 
proposal for the master field can be interpreted as a sort of localisation of the 
functional integral at large- N on a certain microcanonical ensemble specified 
by the set of parameters that describe the master gauge orbit. It is therefore 
natural to classify explicitly gauge orbits. According to the discussion in 
sect. (2) we assume that the relevant gauge orbits are those of hyper-finite 
type, i.e. those which can be obtained as limits of finite dimensional orbits. 
In sect. (4) we will inquire whether and to what extent localisation on this mi- 
crocanonical ensemble is consistent with the loop equations. In this section 
we simply define the microcanonical ensemble. Following our microcanonical 
analogy it is clear that the role of the energy should be played by a set of 
levels of the Hamiltonians that parameterise the gauge orbits. The levels will 
be therefore part of the moduli of the orbits. In addition, since the the gauge 
group acts locally, the moduli will be moduli fields. The representation that 
we are going to find is an identity as it is Eq.(31) and there is a priori no 
guarantee that it will help us to derive an effective action defined over the 
microcanonical ensemble from the loop equations. Indeed it is still necessary 
to functionally integrate over the moduli fields, whose fluctuations a priori 
will contribute loop dependent quantum terms in the loop equations, as it 
happens for the original integration variables, i.e. the gauge connection. In 
the next section we will show that an effective action exists for (planar) closed 
self-avoiding loops, provided we make a further change of variables from the 
levels of the Hamiltonians to the moduli of the associated holomorphic de 
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Rham local system and we allow a central extension in the defining equations 
of the microcanonical ensemble. The existence of the central extension does 
not spoil the hyper-finiteness property in the case of the partially quenched 
Eguchi-Kawai model, but HF is less obvious in the twisted theory since in 
this case gauge orbits have to be infinite dimensional from start. However 
the HF property should be kept also in this case, since the quenched and 
twisted models are equivalent in the large- N limit, though we will not pro- 
vide any direct proof of it. Following our statistical- mechanical analogy we 
want to think the Euclidean functional integral as an integral over the phase 
space rather than the configuration space. DA should be therefore a Liou- 
ville measure associated to some symplectic structure. Since our functional 
integral has a four dimensional nature it turns out that the Liouville measure 
DA is the volume form associated to three different symplectic forms: uj k 
and ujc (ojc is complex, so that its real and imaginary parts provide two real 
symplectic forms): 



In fact the space of four dimensional connections can be thought of as an 
infinite dimensional afline hyper-Kahler manifold. For technical reasons that 
will be evident in sect. (4) we need a reduction of the structure (gauge) group 
from four to two dimensions. This is not restrictive in the large- N limit since 
it is known by Eguchi-Kawai that the theory can be reduced by quenching 
or twisting even to zero dimension keeping the same large- N limit. After the 
partial quenching (twisting) mentioned before, two components of the orig- 
inal connection still transform as a connection under the two dimensional 
structure group while the other two components transform as a Higgs field. 
We notice that the identification of two components of the gauge connection 
with the Higgs field is not canonical, i.e. coordinate independent. However 
in the case of QCD4 we have a preferred choice of coordinates that corre- 
sponds to a flat metric on T 2 x T 2 or better to a conformally flat metric, 
since the classical theory is conformally invariant in four dimensions. Our 
symplectic space can be now identified (non-canonically, i.e. in a coordinate 
dependent way) with the cotangent bundle of two dimensional connections, 
whose coordinates are the pairs (A = A z dz + A 2 dz, ^ = ^ z dz + ^ s dz) . More 




(32) 
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precisely in the partially quenched theory we have the identifications: 



^ 2 = -iD u = -i(d u + iA u )=p u + AZ 
V g = iD u = i(d a - iA n ) = Pu + A°l 



(33) 



while in the partially twisted theory we have: 



\1> = r + A 

= r- + a 



u 



■u 



(34) 



with the other two components of the gauge connection remaining the same 
in both cases. Let us notice that we have used the quenching prescription 
to identify covariant derivatives with the Higgs field in the quenched theory 
and the twisting prescription for the analogous identification in the twisted 
theory. On the cotangent bundle the gauge transformations act by Hamil- 
tonian vector fields with respect to each of the symplectic structures [4, 25]. 
The gauge orbits of the hyper-Kahler reduction are classified by the set of 
levels of the moment maps (the Hamiltonians) for the action of the gauge 
group, one for each of the three symplectic structures: 



with fi° Hermitian and n with values in sl(N). Therefore the partially 
quenched EK model admits a stratification by the following microcanoni- 
cal ensemble: 



where the obvious resolution of identity has been inserted. This equation has 
a meaning for every finite N in the quenched EK model, since the identifica- 
tion in Eq.(33) has a meaning at the level of finite matrices. For the twisted 
model the preceding equations have to be modified because of the non-trivial 
commutator: 




(35) 



Z 




xS(F A - im 2 - ^°)5(d A ij - n)5(d A ifj - n) 



(36) 



u. 



] = 9' 1 ! 



(37) 
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that requires in addition that the gauge orbits be infinite dimensional even 
before taking the large- N limit. In this case our microcanonical ensemble is 
given by: 



modulo gauge transformations, so that a central extension has been intro- 
duced in the Hermitian moment map. From Eq.(38) the microcanonical 
stratification for the partially twisted theory follows: 



The symplectic reduction in Eq. (35,38) and the associated stratification of 
the functional integral are somehow formal, unless the structure of the moduli 
fields is specified. In general, if we require that the symplectic reduction be 
associated to a good moduli problem, i.e. separable, the moduli fields must 
be taken as a linear combination of delta functions, in order that the corre- 
sponding flatness condition defines a (possibly projective) representation of 
the fundamental group of the punctured torus (see appendix B). Despite this 
restriction, these moduli fields are dense in the sense of the distributions in 
the functional integral, that is all that we need to define it. 
In sect. (4) we will see that the existence of punctures, i.e. of a boundary in 
our space-time torus, has a more fundamental meaning to solve loop equa- 
tions than it seems from the definition of the microcanonical ensemble given 
here. We would like to terminate this section writing the resolution of iden- 
tity associated to the microcanonical ensemble in a way that makes manifest 
its four dimensional Euclidean (or Lorentz) covariant nature: 



Z = J DA,D^ U exp(-S YM )S(F^(A) ~ /V " (O^l) (40) 



where the superscript refers to the projection into the anti-selfdual part 
of an antisymmetric tensor. In Eq.(40) we have written the resolution of 




(38) 




(39) 
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identity of the twisted model, but the one of the quenched model can be 
obtained simply setting {0~ v )^ u = and interpreting covariant derivatives 
appropriately. In sect. (5) we will make use of this representation to get an 
explicit Euclidean (or Lorentz) invariant form of the localisation integral. 

4 Solving loop equations by Hitchin systems. 

In this section we consider loop equations for the connection B — A + W 
in the partially quenched or twisted model. The advantage of using B as 
opposed to resides in the fact that planar loops for B generate a complete 
set of gauge invariant observables. The disvantage consists in the fact that 
B is not Hermitian so that ty(x, x; B) is not unitary in general. We will see 
later on under which conditions planar unitary monodromies, i.e. ty(x, x; A), 
can be obtained directly from ty(x, x; B). The loop equations for B read: 

= J DB a exp(-S YM )(Tr(j^^y(x,x;B)) + 

-i ( dy a 5W(z-y)Tr(*(x,y;B))Tr(*(y,x;B))) (41) 

JC(x,x) 

with a = 1, 2, so that in these variables they are as difficult to solve as for the 
original four dimensional connection A^. The quantum contribution in gen- 
eral is loop dependent and this dependence is an obstruction to the existence 
of a conventional (hyper-finite?) effective action (we take the point of view 
that an effective action defined directly in loop space [48, 49], from which 
the loop equations may be derived as critical equations, probably would not 
be hyper-finite). 

After implementing in the functional integral the resolution of identity by 
means of the gauge orbits of the microcanonical ensemble of sect. (3), we 
change variables in the loop equations from Hitchin bundles to the corre- 
sponding holomorphic de Rham local system (see appendix B). In [59] an 
analogous change of variables, from the Hitchin system to the action-angle 
variables of the Hitchin fibration of the corresponding Dolbeault moduli space 
(see appendix B), was proposed on the basis of energy /entropy considera- 
tions. 

Our basic philosophy consists in changing variables in such a way that the 
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quantum contribution may be computed explicitly by evaluating a residue 
that depends only on the topological class of the loop. It is well known by the 
Cauchy formula that the line integral along a closed loop of a holomorphic 
function times the Cauchy kernel with pole not lying on the loop depends 
only on the winding number of the loop around the pole and on the value of 
the holomorphic function at the pole of the Cauchy kernel: 



where Indc(z) is the winding number of C around z. The existence would 
become possible of an effective action, for loops in a certain topological class 
determined by the winding number, if only the Cauchy kernel could be gen- 
erated by functionally differentiating the monodromy of the connection in 
the loop equations with respect to the integration variables in the functional 
integral. We have seen in sect. (2) that functional differentiation with respect 
to the gauge connection produces a delta function in the quantum contribu- 
tion to the loop equations. The change of variables to the de Rham systems 
is thought in such a way that functional differentiation of the monodromy 
with respect to the moduli fields of the holomorphic de Rham local systems 
produces precisely the Cauchy kernel in the loop equations. However, even 
if we reduce the computation of the quantum contribution to the evalua- 
tion of a line integral of a Cauchy kernel, there are two obstructions for the 
residue theorem to apply. The first one is that the monodromy is not a holo- 
morphic function of its endpoints. The second one is that gauge invariance 
requires that the functional derivative with respect to the integration vari- 
able be taken at a point of the loop, for the result to be non-vanishing, i.e 
at a singularity point of the Cauchy kernel. We will see below how each of 
these obstructions is the cure for the other one. After introducing the non- 
Hermitian connection B = A + i^fr = b + b the microcanonical representation 
of sect. (3) holds in the form: 



where we have re-casted the microcanonical resolution of identity in a dif- 
ferent, equivalent form (see appendix B), for later convenience, in order to 
compute the effective action (sect. (5)). We now change variables from the 




(42) 



Z = / 5(F B - ii)5(d* A y - v) exp(-S EK )DBDfiD 



v 



(43) 
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Hitchin bundles to the holomorphic de Rham local systems (see appendix 
B), whose moduli are described by the field //, obtained from the equation: 

F B - ll = (44) 

by means of a complexified gauge transformation G(x; B) that puts the con- 
nection B = b + b in the holomorphic gauge b = 0: 

Bb - i// = (45) 

where li = fjP + (n — n) and // = GliG^ 1 . The partition function in terms of 
the moduli fields of the de Rham local system becomes: 

Z = J 5(F B - Li)5{d* A m - v) exp(-S EK )DB ^ ^ Dji'Dp! (46) 

The integral over B can now be performed and the resulting functional de- 
terminants, together with the Jacobian of the change of variables to the de 
Rham local system, absorbed into the definition of the effective action T. T 
will be written more explicitly in the next section in terms of these functional 
determinants. Let us notice that for computing functional determinants a 
choice of metric is necessary The natural metric for us is the translational 
invariant metric on the torus exp(— p)dzdz with the conformal factor con- 
stant. However the classical action, being conformally invariant, does not 
depend on the choice of p. Yet we will see below that our regularization of 
the quantum contribution to the loop equations will require to transform the 
metric by a conformal diffeomorphism. The partition function is now: 

Z = J exp(-r)Dn'Dji' (47) 

To realize our aim of getting the Cauchy kernel by functional differentiation 
with respect to the moduli fields of the de Rham local system it is convenient 
to study the loop equations for the Wilson loop constructed by the connec- 
tion b, thought as a functional of B corresponding to gauge transforming 
B into the gauge b = 0. Such a gauge transformation belongs to the com- 
plexification of the gauge group and it is rather a change of variables than a 
proper gauge transformation. However, because of the property of the trace, 
for closed loops, it preserves the trace of the monodromy. This will allow us, 
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later in the course of our computations, to transform the loop equations thus 
obtained into equations for the monodromy of B. In our derivation of the 
loop equations, a crucial role is played by the condition that the expectation 
value of an open loop vanishes. For open loops constructed by means of the 
connection B this property follows directly from the local gauge invariance 
of the partially quenched or twisted EK model. However, the starting point 
in our case is the monodromy of the connection b, whose expectation value 
is also required to vanish for open loops. Yet, this condition does not follow 
from gauge invariance, since the change of variables to the de Rham local 
system (b = 0) does not allow any residual local compact gauge freedom (un- 
less perhaps we analytically continue to Minkowski space-time, see below). 
We will present here two slightly different derivations of the effective action 
each based on a different way of achieving the vanishing of the expectation 
value of open b loops. 

The most direct procedure consists in observing that in the partially quenched 
theory the vanishing of the expectation value of open b loops follows from 
the fact that the partial quenching is thought precisely to keep unbroken the 
R 2 symmetry B — > B + i(p u dz + pudz) and that in the gauge 6 = the R 2 
symmetry is broken down to b — > b + ip u dz. If this residual symmetry is 
unbroken, then the expectation value of the trace of open b loops vanishes 
because of the equation: 



N-^TrViz!, z 2 ; b)) = exp(p tt (z 1 - z 2 ) z )N~ 1 (Tr*( Zl , z 2 ; b)) (48) 



provided the normalised trace stays uniformly bounded in the large- N limit. 
To be precise, this holds only for loops that are open once projected on the 
z direction, because the R 2 symmetry is partially broken to b — > b + ip u dz. 
However this weaker property is all that is needed in our loop equations. In 
the other proof of the loop equations that we propose below, the assumption 
of uniform boundness is not necessary and the vanishing of the expectation 
values of open b loops holds without further specifications. We may thus 
derive our loop equations in the partially quenched theory: 
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+ / dy z d-\w-y)Tr{\ a m{x,y ] b)\ a m{y,x ] b))) = 



JC(x,x) 



= fDfjL'Dji'exp(-r)(Tr(— — *(x,x;b)) + 
J dfx'yw) 

+ / dy z d- 1 (w-y)Tr(^(x,y;b))Tr(^(y,x;b)) + 



JC(x,x) 



-±-Tr{V{x,y;b)y(y,x; &))) 



(49) 



that in the large- iV limit reduce to: 



0= f Dfi'Dp!exp(-r)(Tr(— --$>(x,x;b)) + 
J o/i(w) 

+ [ dy^'^w - y)Tr(^(x,y;b))Tr(^(y,x;b))) 



JC{x,x) 



(50) 



The most interesting case is when w lies on the loop C, because the loop 
equations can be transformed easily into an equation for B. However in 
this case the contour integration includes the pole of the Cauchy kernel. 
We need therefore a gauge invariant regularization. Before doing so it is 
convenient to make a conformal diffeomorphism to put the contour integral 
in standard form. We wish to do that because, when we regularize the theory, 
we require a loop-independent universal regularization. Hence it is natural 
to use the uniformization theory of Riemann surfaces with a boundary. It 
is not restrictive to take w as the base point of the loop. For every simple 
loop (i.e. non-overlapping, self-intersecting only at the base point) we make 
a conformal transformation that maps each petal into a disk and then each 
disk into the upper half-plane. Each petal of the loop is then mapped into 
a copy of the upper half plane and the boundary identified with the x-axis. 
Let t be the local coordinate on the upper half-plane corresponding to this 
conformal transformation. Obviously t depends on C and the microcanonical 
resolution of the identity in the functional integral, once written in term of 
t, depends on C too. The loop equations in the variables associated to the 
coordinate t look like the ones in the original global coordinate on the torus 
z. Therefore we get: 







J D^ f L>^-exp(-r)(Tr( 



5T 




(51) 
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with the important difference that the contour integral is now restricted to 
the real axis and that the metric needed to define the functional determinants 
in T is the conformally transformed metric exp(—p)dtdt = exp(— p)\^\ 2 dzdz. 
In d = 4 it is expected that this conformal transformation will affect the ef- 
fective action only through its dependence on the conformal anomaly [60] . A 
natural regularization of the quantum contribution to the loop equations is 
a ie regularization of the Cauchy kernel. This is the most natural among any 
other possible regularization in Euclidean space, since it is in fact compatible 
with the analytic continuation to Minkowski space-time, that manifestly pre- 
serves gauge invariance. The idea that Feynmann diagrams in two dimensions 
could be regularized by analytical continuation to Minkowski space-time in 
certain cases is contained in [57]. A master field living in the Cuntz algebra 
in Minkowski space-time has been constructed in general, and in particular 
for Yang-Mills theories, in [61]. The result of the ie regularization is the sum 
of two distributions, the principal part plus a one dimensional delta function: 

d -1 ^ - y x + ie) = {2ti)- 1 (P(w x - y x y l - ind(w x - y x )) (52) 

The loop equations thus regularized look like: 

0= [ Df i'Dfi>ex P (-r)(Tri-^— *(x,x;b)) + 
J dn'(w) 

+ / dy x {2i:)~ 1 (P(w x - y x )~ l - m5(w x - y x )) x 

JC(x,x) 

xTr(^(x,y;b))Tr(^(y,x;b))) (53) 

where we are omitting the t under-script. Being supported on open loops the 
principal part does not contribute and the loop equations reduce to: 

r 5T 
= j Di2'Dp'exp(-T)(Tr(j-^-^(x,x-,b)) + 

- I dy x l -5{w x - y x )Tr(V{x, y; b))Tr{^(y, x; b))) (54) 

JC(x,x) Z 

Taking w = x and using the gauge transformation properties of the b mon- 
odromy and of /u(x)', the preceding equation can be rewritten in terms of the 
connection B and the moduli fields /x: 

r 5T 
= j Dp'Dp'exp(-T)(Tr(j-^^(x,x-,B)) + 
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(55) 



where we have used the condition that the trace of open loops vanishes to 
substitute the b monodromy with the B monodromy. This is our final form 
of the regularized Euclidean loop equations. Let us notice that there is a 
sign ambiguity in the quantum contribution that depends on the choice of 
the sign of e in the regulator or as to whether the conformal transformation is 
orientation preserving or orientation inverting. The occurrence of the factor 
of | seems to be related to the fact that the pole of the Cauchy kernel lies 
along the contour integral. If the contour integral were to include the pole 
inside the loop and the residue theorem to apply naively, the residue would 
be 1, while it would be if the pole were outside the loop. Both choices 
however would violate gauge invariance. 

For general self-intersection the Euclidean loop equations reduce to (assuming 
the conformal map to be orientation preserving): 



Because of the product of traces, in general the quantum contribution does 
not have the same operator structure as does the classical term and it is loop 
dependent. Therefore, in general, an effective action need not to exist as 
a functional over the microcanonical ensemble constructed by the inductive 
limit of finite dimensional orbits. However, for self-avoiding loops, one of 
the traces in the quantum contribution is just Tr(l) = N and there is only 
one term as for the classical contribution. It is now clear the role played 
by the absence of self-intersection for these variables associated to the de 
Rham system of the microcanonical ensemble. It is only in this case that the 
quantum contribution is proportional to the monodromy as is the classical 
term. In general we get a loop dependent factorized trace of monodromies, 
that cannot be re-absorbed into a redefinition of T. 

The result is the following Euclidean master equation for self-avoiding loops: 



r dl 

= J D^D/2'exp(-T)(Tr(j-^ty(x,x-,B)) + 

- % - J2 Ind c (y)Tr(V(x,y;B))Tr(y(y,x-B))) 



(v;v=x) 



(56) 



J exp(-r)(7V(( 



sr 



l -Ind c (x)Tr(l))*(x,x;B))) (57) 



= 



Sfi{x) 2 
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that reduces to: 



5T 



l -Ind c {x)Tr{\) 



(58) 







Sfi(x) 2 



The corresponding effective action is: 



Y eff = Y-e{C) l -jTr^. 



x 



(59) 



with e(C) = +1 or — f. The two possible signs are due to the two possible 
orientations of the loop C. In addition there is a sign ambiguity due to the 
conventional choices that we made: the choice of regularization, i.e ie versus 
— ie and the convention that the conformal map is orientation preserving. 
Remarkably the sign ambiguity in the Euclidean loop equations corresponds 
in fact to physical equivalence. Indeed the quantum contribution is a ^-term 
with 8 — \. Now, since / Tr/j,d 2 x represents the first Chern class, 8 is defined 
modulo integers, so that for 8 — | a change of sign does not change 8 modulo 
1. Had we regularized by analytical continuation to Minkowski space-time 
we would have obtained the same imaginary quantum contribution to the 
effective action, while the action T would have been Wick rotated by a factor 
of i. In any case we see that the quantum contribution for the de Rham 
system has produced a central extension in the master equation. This means 
in general, that even if we start with a microcanonical ensemble associated 
to sl(N) bundles, central directions in the moduli fields have to be included, 
otherwise the master equation cannot be satisfied. Hence, in the partially 
quenched model we have to extend our microcanonical ensemble from sl(N) 
to gl(N) and to allow a central extension in the Hermitian moment map as 
follows: 



After that the entire chain of derivations of this section can be consistently 
performed to arrive at the master equation. Yet, the occurrence of the 6*-term 
in the effective action is not completely harmless, since it affects the topo- 
logical sectors of the U(N) theory (see for example [62]). In fact, we will see 
in sect. (5) that the central term in F e ff is cancelled by a counter-term, that 



Fa - i^ 2 = fiP + cl 

dAip = n 



(60) 
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arises as a -u(l)-anomaly in the Jacobian of the "chiral" change of variables 
from the Hitchin to the holomorphic de Rham systems, after a finite con- 
formal rescaling of the metric used to compute the functional determinants. 
The "holographic" nature of our change of variables will be then completely 
manifest, in the sense that we are cancelling the quantum boundary contri- 
bution against one from the bulk, to obtain an effective action defined purely 
on the bulk. 

Summarising, given a sequence of solutions of the master equation for any 
finite N, we have a sequence of solutions of the loop equations for finite N, 
up to terms vanishing in the large- N limit, corresponding to the omitted 
connected contribution. Therefore for iV large but finite we have an ap- 
proximation of the solution of the N = oo bare loop equations for planar 
self-avoiding loops with accuracy iV -1 . The algebra generated by the mon- 
odromy of the master connection of the partially quenched EK model for 
any finite iV is a matrix algebra, since the microcanonical ensemble defines a 
finite dimensional representation (of a central extension) of the fundamental 
group of a punctured torus. Obviously renormalization is necessary, since we 
have written the master equation for the bare master field, but in principle 
nothing forbids to perform it. Assuming convergence of the renormalised 
master monodromy as N goes to infinity, it follows the hyper-finiteness of 
the master algebra for planar closed self-avoiding B loops. 
It is also clear that the twisted model leads to the same master equation, 
with the only difference that now the central extension is present from start 
and it is associated to an infinite dimensional projective representation of the 
fundamental group. 

To show hyper-finiteness also in the twisted case we must associate to this 
infinite dimensional representation a system of approximating finite dimen- 
sional representations. In the twisted case it is not restrictive to start with a 
microcanonical ensemble that, though is infinite dimensional because of the 
central extension due to the non-trivial commutator: 

[r u ,r u ] = 9- 1 i (6i) 

has finite rank residues at the punctures. This implies that the monodromies 
at the punctures are a direct sum of finite dimensional matrices and a multi- 
ple of identity. The representation obtained projecting away the multiple of 
identity part is the desired approximating system (finite dimensional approxi- 
mations of infinite dimensional representations of the fundamental group have 
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been previously considered in [63] in the framework of non-commutative Rie- 
mann surfaces and in the context of [59]). To be explicit the master algebra 
is the algebra generated by an inductive sequence of finite dimensional group 
algebras of representations of the fundamental group, where the inductive 
embedding is with respect to a partial order on the growing divisor, D, and 
on the increasing rank, N, of the local system. Thus we get in the quenched 
case: 

*(a)*(6)*(a)- 1 *(6)- 1 JJ V(D p ) = 1 (62) 

v 

where the product is over the monodromies around the punctures and a 
and b are the cycles of the torus T 2 . We assume translational invariance, 
that means that the monodromy around the punctures are constant up to a 
possible conjugacy transformation by means of the operator of translations 
of the quenched theory. In the twisted case the representation is projective 
and the constant monodromies may be conjugated either by the translation 
operator or by another twist operator: 

*(o)*(6)*(o) _1 *(6)- 1 n*Pp) = exp(*cV (2) )l (63) 

v 

where is the area of our space-time torus. The next version of the 
derivation of the effective action from the loop equations is based on the fact 
that for U (N) bundles, that we have in any case to introduce because of the 
occurrence of the central extension in the master equation, we can realize 
explicitly the R 2 symmetry as a global symmetry, even for the de Rham 
systems. In fact what really matters to generate the Cauchy kernel in the 
loop equations is the relation: 

Bb -in' = (64) 

that follows from Eq.(35) in the gauge 6 = 0. Eq.(64) holds under slightly 
weaker conditions than 6 = 0, that allow us to keep explicitly a residual R 2 
global symmetry. For the Lie algebra u(N) we decompose the connections 
B and b into their semi-simple and central parts: B — B s + B c , b = b s + b c . 
Our gauge fixing is now: 

db c = 

I s = (65) 
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This gauge fixing admits the residual global symmetry b — > b + ip z dz, b ^b — 
ip z dz maintaining at the same time the validity of Eq.(64). If this symmetry 
is unbroken, than the expectation value of the trace of open b loops vanishes 
because of the equation: 



(ty(x,y;b c )) = exp(ip(x - y))(^(x,y;b c )) (66) 



We get: 



0= f (fpDfi'Dfi'Tr - - . (exp(-T)V(x,x;b)) (67) 

J dfJ,{W) 

where we have inserted the integration over the zero mode d 2 p. In the same 
vein as before we obtain: 

= / d 2 pD^'Dji' exp(-r)(Tr(-^—^(x,x;b)) + 
J dfi'{w) 



+ / dy z d-\w-y)Tr(\ a m(x,y-b)\ a m(y,x ] b))) = 

JC(x,x) 

= [ d 2 pDfi'Dji'exp(-T)(Tr(-^—y(x,x;b)) + 
+ / dy z B~\w - y)Tr(*(x, y; b))Tr(*(y, x; b))) 

JC(x,x) 



(68) 



Then the master equation follows as before. Several observations are in order. 
An important feature of our regularized master equation is that it depends 
implicitly on the conformal diffeomorphism that maps the interior of the 
loop into the upper half plane. This is a direct consequence of our choice 
of regularising the quantum contribution to loop equations in a universal 
way. Indeed the local field for which the loop equations can be reduced to 
the master equation is the field fi t t that transforms as a (1,1) form under 
conformal diffeomorphism. Instead /j zz is the field that occurs in the classical 
action, where (z, z) are the global coordinates on the torus. The existence of 
this conformal mapping plays a key role when we introduce a lattice divisor, 
D, to get an explicit realization of our Hitchin system, actually the only one 
possible if we require a Hausdorff moduli space. In this case \x is a collective 
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field given, according to sect. (3) and appendix B, by: 



H = ^4S^(x-x p ) 



p 



n = ripd^ (x — x p ) 



p 



n = n p 5^ (x — x p ) 



(69) 



p 



Thus ^ is a "lattice (or adelic) field of residues". As a consequence we need a 
more sophisticated version of the uniformization map. In this case we ought 
to map conformally the interior of the loop minus the internal punctures of 
the lattice divisor into a punctured disk or into a fundamental region in the 
upper half-plane. In this case the image of the punctures in the upper half 
plane are cusps touching the real axis. The loop equations in terms of the 
"lattice field of residues" are thus: 



The cusp singularity at the point x p introduces one more sign ambiguity for 
the central term in the effective action. In the lattice field version of our loop 
equations, represented by Eq.(70), we must take by necessity the base-point 
of the loop at a puncture, otherwise we get a trivial result. This may seem 
unnatural from the point of view of the uniformization map, but somehow 
surprisingly this is precisely the presentation that occurs in Penner [64] defi- 
nition of the Teichmuller space of bordered surfaces, whose moduli involve an 
arbitrary choice of a point on the boundary, that is then projected to infinity 
to become a cusp on the boundary of the upper-half plane. In appendix C 
we will see how the existence of the conformal diffeomorphism involved in 
the derivation of our master equation for the field of residues seems to be 
crucial to get the correct area law for the Wilson loop in two dimensions by 
the Hitchin systems. 

In a sense we may consider the effective action as the image of a "holo- 
graphic" map from the boundary, i.e. the loops, to the bulk, i.e. the torus 
T 2 , the space-time of the quenched theory. From this point of view the HF 
property of the master algebra is related to how much entropy is encoded 




(70) 



29 



into the boundary. For self-avoiding loops we have a HF effective action, 
while for arbitrary self-intersection of the boundary we may need a non-HF 
algebra. 

The preceding derivation of the master equation from the loop equations lit- 
erally allows us to employ the master field B, that is the solution of Eq.(70), 
only to compute the corresponding Wilson loop limjv N~ l Tr^(x, x; B). How- 
ever the quantity that has most physical interest is rather lim^r N~ 1 Tr^{x, x; A), 
i.e. the trace of the unitary monodromy. This quantity can be computed via a 
limiting procedure taking advantage of the existence of the following k action 
(ip — > Kip) on the Hitchin bundles [5] (see appendix B). After some obvious 
rescaling in the defining equations of our microcanonical ensemble we get a 
k dependent effective action appropriate for the study of the monodromy of 
the operator B K = A + i(nip + Rip) = A K + i^ K : 



After that the monodromy of A is recovered in the limit k — > 0. This poses 
the interesting problem as to whether there exist solutions for A and ^ that 
are k independent, as it would be required by the spirit if not the letter 
of the master field philosophy, and more generally Euclidean (or Lorentz) 
covariant, in the sense that the Euclidean or Lorentz group are embedded 
into the (global) gauge group [10]. As far as ^-independence is concerned 
it is clear that any ^-independent solution must coincide with a fixed point 
of the k action as k — > on the infinite dimensional Hitchin bundles in- 
volved. The study of the limit k — > would also furnish the Wilson loop 
for the A monodromy. We are therefore interested to study the k action 
on infinite dimensional Hitchin systems, a case on which we have not much 
information. However in finite dimensions such fixed points exist and form 
a Lagrangian subspace of the moduli of Hitchin bundles that coincides with 
those bundles that arise as variations of Hodge structures (see appendix B) 
[65, 66, 67]. For such bundles the B monodromy coincides with the A mon- 
odromy [65, 66, 67, 68] (see appendix B). As far as Euclidean or Lorentz 
covariance is concerned, bundles that arise as variations of Hodge structures 



Z = J 5(F BK -fi- ^ 1 |k| 2 1)5(^ k * k -u)x 
x exp(-S TW )DB ^; U I ds = 




(71) 
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are associated to embeddings of sl(2, R) into sl(N, C) via the Schmid nilpo- 
tent orbit theorem [65, 66, 67]. It is natural to speculate, taking into account 
also the peculiar interplay between the Euclidean and Minkowski version of 
the master equation, that such embeddings may be related to embeddings of 
the Lorentz group into the gauge group. In addition we observe that a neces- 
sary condition for Euclidean or Lorentz covariance is given by the vanishing 
in the large- N limit of the central extension in the microcanonical ensemble 
for the master field. Indeed any central element can only transform as a 
scalar under Euclidean or Lorentz transformations, if these are to be embed- 
ded into the gauge group, while this contradicts the nature of the defining 
equations of the microcanonical ensemble, that transform as a anti-selfdual 
tensor in our four dimensional interpretation. We have already seen that the 
effective action develops a central term as result of the quantum contribu- 
tion. However this term must be an artifact due to the "chirality" of our 
change of variables. Therefore it must be cancelled by the contribution of 
an "anomalous" Jacobian of the change of variables in the bulk. In the next 
section we will see that T contains a counter-term that cancels the quantum 
contribution provided a suitable finite rescaling of the conformal factor in the 
metric is made (indeed the " anomalous" Jacobian depends on the conformal 
factor of the metric, once the complex structure is fixed, while the quantum 
contribution to the effective action, being a residue, does not.) 



5 The effective action of large- N QCD4 and 
of its Af = 1 extension 

We are now ready to look at the effective action for Hitchin systems in some 
detail. It easy to see that the effective action, as it would be computed from 
the loop equations without gauge fixing, is the sum of the classical (quenched 
or twisted) Yang-Mills action plus the logarithm of a number of determinants 
whose origin is the following one. 

First of all there are the determinants that arise from the localisation integral 
on the microcanonical ensemble, i.e. from integrating the delta functionals 
in Eq. (36,39) over the space of connections: 

Det'{d B A)~ 1 Det'{Ud* B )~ 1 Det'ip*uj) (72) 
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The first two determinants arise from integrating over the tangent space to 
the Hitchin equations. The third one is the pull back, on the moduli that 
are not local moduli, of the symplectic form on the space of connections. 
According to the formulas on the dimension of the moduli space in appendix 
B, there are not such non-local moduli on a torus or on a sphere, once the 
local moduli fields are fixed, while there are on Riemann surfaces of higher 
genus. Therefore for the quenched or twisted theory on a torus the third 
factor is absent. The ' suffix refers to the need of projecting away from 
the determinants the zero modes due to gauge invariance, since gauge-fixing 
is not implied in the loop equations, though it may be understood if we 
like to. Finally there is the determinant that arises as the Jacobian of the 
change of variables from the local moduli fields of the Hitchin system to the 
corresponding holomorphic de Rham local system: 

^#r4 (73) 

At formal level and without an explicit gauge fixing, the effective action, T, 
reads: 

T = S YM - logDet'(d B A) - logDet\Ud* B ) + logDet' f*^' ^ (74) 

where Sym is the Yang-Mills action of the reduced quenched or twisted the- 
ory. The effective action can be given a non-formal meaning introducing a 
metric, to make sense of the formal determinants, and making a choice of 
gauge by the Faddeev-Popov procedure, to avoid the zero modes. A conve- 
nient choice of gauge is [69]: 

d* A 8A + iadlS^ = (75) 

that coincides with a background Euclidean or Lorentz covariant gauge (after 
analytical continuation to Minkowski space-time). This gauge fixing condi- 
tion can be written in terms of B as [69]: 

U{d* B 5B) = (76) 

After introducing this gauge-fixing and the corresponding Faddeev-Popov 
determinant, the localisation integral over the microcanonical ensemble gives 
us: 

Det(-A 1 (S))-5 J D e t(_sRA B )'^ J Det(-SA B )'5 J Det(_KA B )- 1 (77) 
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where is the scalar Laplacian and Ai(B) is the operator: 

Ai(S) = d* B d B A (78) 

acting on one-forms and restricted to d* B 5B = 0. The operator —Ai(B) 
arises from interpreting / DB5(Fb — n) as: 

lim / DBex V (--Tr((F B - p)(F B -//)*)) (79) 

with the constraint d* B 5B = 0. The last determinant in Eq.(77) is the Faddev- 
Popov determinant. The last three factors in Eq.(77) cancel together. The 
effective action is thus: 

r = S YM - hogDeti-A^B)) + logDet (80) 

The first determinant, that we obtained from the microcanonical localisation, 
is essentially a one-loop contribution in a background field. We should notice 
that it contributes two Hermitian polarisations. In fact the connection B is 
non-Hermitian and thus it is determined by four real polarisations, but two 
of them are killed by the conditions: 

fft(cT B SB) = 

%(d B 5B) = (81) 

that come from the gauge-fixing and from the harmonic constraint [69] in the 
functional integral. However this determinant can be more conveniently de- 
fined and computed by means of a manifestly Euclidean (or Lorentz) invariant 
metric and of the representation of the microcanonical ensemble reported at 
the end of sect. (3). Though we use now a manifestly Euclidean (or Lorentz) 
invariant notation, there is no problem to interpret the following formulas 
appropriately for the quenched or for the twisted models. The localisation 
integral in Eq.(40) is defined formally as: 

J DA,6(F^ U - fi- u ) = Det'-\p-d A /\) = Det'-^-A,) (82) 

where P~ is the projector onto the anti-selfdual part of the field strength. 
After inserting the gauge-fixing condition and the corresponding Faddeev- 
Popov determinant it can be defined non-formally as: 

i t 1 

Det~3(-Ai) = lim / DcDA^ exp( — Tr(c 2 )) x 
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x eM A Tr ((F- - ^ v f)5{D,8A, - c)A FP (83) 

The result for Det~^{— Ai) is thus: 

Det~3(-Ai) = Det~2(-AA<5^1 + iad F -)A FP (84) 
where is the anti-selfdual part of the field strength, given by: 

F» v = F »v ~ F tv (85) 

with: 

F lv = -^vvapFaP (86) 

It is a fundamental problem to study the renormalization group properties 
of the master equation for our effective action, with special regard to the 
asymptotic freedom [70, 71]. In appendix D we have computed the divergent 
part of the localisation determinant Det~^{— Ai) considered as a functional 
of the anti-selfdual F~ u or self-dual field strength F+ u : 

1 11 TV A 2 r 1 

--(ojBetf-AO = --— to 9 (— ) J d'x-Tr(F-) 2 + 

4<l^ ( f>/ <fttTr < F -> 2 < 87 > 

We should notice that the logarithmic divergence associated to Tr(F+ y ) 2 
really occurs in the effective action only if the effective action is regarded as 
a local functional of F~ v and F+ v . Yet in the master equation we should rather 
consider the effective action as a functional of the moduli fields. Now the 
components of the tensor F~ v are certainly a local functional of the moduli 
fields, since the moduli fields are a linear combination of them, while F+, is 
not, since: 

p-{d A 8A) = SfT (88) 
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Because of this non-local character it is not obvious that, once the non-local 
expression in terms of the moduli fields is inserted in the definition of the 
localisation determinant, it produces a logarithmic divergence. In any case 
the master equation reads: 

?¥{ = (89) 

The left hand side contains, by the definition of the moduli fields //, a linear 
combination of: 

T^fr 00) 

Using Eq.(88) we get Jd 4 xTr(F^) 2 = / d A x\Tr{F~) 2 at lowest order in 
powers of n~ v , assuming Euclidean (or Lorentz) invariance (of the quenched 
theory). It follows that, at first order in the expansion in powers of the 
moduli fields, which by the way coincides at this order with the usual per- 
turbative expansion, the contribution of the self-dual part of the effective 
action vanishes. The first coefficient of the beta function is thus exactly re- 
produced at this order by the localisation determinant (for the moment we 
do not include the contribution of logDet -§^r^ whose computation we leave 
for the future): 

J^j^~Tr ( F^ - %J^-M$)S ^-Jr(F^) = (91, 

The preceding construction of the effective action of QCD4 extends easily 
to the super-symmetric case without scalars. In fact, if only fermions are 
present in addition to the gauge fields, they can be integrated out, the result 
being a Pfafnan whose logarithm should be added to the effective action. 
Thus in the Af = 1 case, we get for the effective action: 

T = S YM -\logDet(-A 1 (B)) + 

+ l °9Det ^; U J + logPf(aM (92) 

We must now give an interpretation of the conformal diffeomorphism, due 
to the uniformization map, compatible with the four dimensional interpre- 
tation of the quenching or twisting procedure. Thus the dependence on the 
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conformal factor in the first determinant of Eq.(80) will be only through 
the conformal anomaly, since it is the determinant of a Laplacian acting on 
one-forms in four dimensions (though the Laplacian looks two dimensional, 
we have seen that the missing two dimensional integration is hidden into 
the definition of the trace and in the definition of \1/ in the twisted theory or 
into the dimensionful rescaling of the coupling constant and in the quenching 
prescription for \I/ in the quenched case). The second determinant in Eq.(80) 
instead depends explicitly on the conformal factor of the metric since it is 
the determinant of an operator acting on scalars. One of the features that 
distinguishes the present approach from the microcanonical localisation, is 
the contribution from the Jacobian of the change of variables from Hitchin 
bundles to the de Rham systems: 

D(fj,,u) _ x DQSjjT) _ D i ( % %\ m) 

where to get the zero in Eq.(93) we used the relations: 

d* B 5B = ibv 
d B A5B = 5fi (94) 

This determinant is in fact an " anomalous" Jacobian for the action of the two 
dimensional complexified gauge group. This is reminiscent of the fact that 
the analogous Jacobian that occurs in the collective field method applied to 
loop equations arises as an anomaly in the algebra of loop substitutions [72]. 
More explicitly we can extract from logDet ®^; u ^ a term linear in Trfi as 
follows: 

DeA = Der*M) = 

= Det~\^^-) = D ^ + G - 1 ^) (95) 

We now use the formula: 

Det(l + A) = exp Trlog(l + A) (96) 

Since A in our case is a commutator, in finite dimension the contribution 
linear in A would vanish. However this is not anymore true in infinite di- 
mensions, that is our case, since the trace includes integration over space- 
time. We will apply the general theory of regularized determinants to the 
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case that A is a pseudo-differential operator. According to the general the- 
ory [73, 74], the definition of regularized determinant employs a generalised 
noncyclic trace TR reg . The definition of TR reg involves a choice of a positive 
pseudo-differential operator Q of positive order q. 

TR re JA) = MmTr(exp(-eQ)A) (97) 

e— >0 

or equivalently (the two definitions are related by a Mellin transform): 

TR reg (A) = lim Tr(Q- z A) (98) 

TR reg is a non-cyclic trace, so that, in general, the trace of a commutator 
does not vanish. The renormalised trace is instead: 

TR ren {A) = lim Tr(Q' z A - ^-res(A)) (99) 

where res(A) is the Wodzicki residue of a pseudo-differential operator. It is 
defined as: 

res ( A ) = / dx I tra- d (x,p) (100) 

{27T) a JM J\p\=l 

where a-d(x,p) is the homogeneous term of order —d in the momenta, in 
the asymptotic expansion of the matrix valued symbol of A. The Wodzicki 
residue corresponds to extract the coefficient of the logarithmically divergent 
term in the symbol expansion of a pseudo-differential operator. Since usually 
anomalies arise from finite integrals, for anomalies the regularized and renor- 
malised traces coincide. An important formula to compute the regularized 
trace is: 

TR reg [A, B] = ~res([log\Q\,A]B) (101) 

We are now ready to discuss the contribution of the chiral Jacobian Det(jjjj). 
We anticipated that in the effective action the central term that comes from 
the residue in the quantum contribution to loop equations must be cancelled 
by an analogous term from the chiral "anomalous" Jacobian, since QCD is 
non-chiral and there is no central term in the defining classical theory. This 
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anomalous term in the chiral Jacobian, that is the trace of the term linear in 



I DC 

-resdloglQlG- 1 ^) 



(102) 



defines a cocycle, that coincides with the Radul cocycle. Now the commu- 
tator anomaly for point-wise multiplication vanishes in every dimension but 
one, where coincides with the (twisted) Radul cocycle. However our com- 
mutator involves an operator, G -1 -^-, whose symbol is not a local operator, 
since the functional derivative of the gauge transformation with respect to 
fj, involves the inverse Laplacian. Thus in our case the commutator anomaly 
does not vanish in dimension two, if we take as Q the two dimensional Lapla- 
cian, a choice that seems appropriate in the quenched theory. Thus we may 
compute the u(l) part of our chiral anomaly as follows (for the metric dzdz): 



From Eq.(103) we see that there exists a finite mismatch of a factor of n 
between the commutator anomaly and the quantum central term in the ef- 
fective action. However the commutator anomaly depends on the choice of 
the metric, in particular of the conformal factor of the metric, while the met- 
ric and a fortiori its conformal factor do not affect the residue in the quantum 
term. Therefore we are able to cancel the central term in the effective action 
in the quenched theory by a conformal rescaling of the metric of a factor 
of 7r. It would be interesting to perform the analogous computation in the 
twisted theory. In four dimensions a finite conformal rescaling of the metric 
changes the effective action only by finite terms due to the four dimensional 
conformal anomaly. The same happens in two dimensions through a differ- 
ent mechanism though: since the classical action depends on the conformal 
scale of the metric, a conformal rescaling is equivalent to a rescaling of the 
dimensionful coupling constant (see appendix C). To be precise the sign of 
the commutator anomaly in the Euclidean space-time is fixed, while it is not 




(103) 



38 



so for the sign of the quantum contribution, so that the cancellation between 
them can actually occur only modulo integers, but in fact this is all that is 
needed, since the 6>-term in the effective action is defined modulo integers 
too. 



6 Loop equations and the geometric Lang- 
lands correspondence. 

We have seen in the previous section that loop equations admit as solutions 
inductive limits of (projective) representations of the fundamental group of 
a punctured torus. In this section we will use the geometric Langlands cor- 
respondence in the sense of Beilinson and Drinfeld [26, 27] (and references 
therein) to interpret the master equation by operator solutions of the loop 
equations obtained by the quantisation of the Hitchin systems. As we will 
see the construction that we are referring to is in a certain sense dual to the 
one of sect. (4). Because of quantisation, the HF of the ambient algebra of the 
master field is not preserved a priori, while it is kept quantum integrability. 
We have seen that the loop equations: 

0= [ D/i'Dfi'exp(-T)(Tr(-^—^(x,x;B)) + 
J dfj,{w) 

+ f dy z d-\w - y)Tr(*(x, y; B))Tr(*(y, x; B))) (104) 

JC(x,x) 

for self-avoiding loops are implied by the master equation: 

= f$ (105) 

However there is a different way of satisfying the loop equations. We may 
think that the functional integration in the loop equations defines a gener- 
alised trace (Q 1 1 Q) on the field algebra. Quantisation in presence of con- 
straints, that in our case are Hitchin self-duality equations on a punctured 
torus, is usually performed in two ways. Either we first quantise the space 
of connections [35, 36] and then we impose the constraints on the physical 
states of the theory or we first solve the constraints and then we quantise the 
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reduced theory. In our case we must follow the second alternative because 
of the localisation over the microcanonical ensemble in the functional inte- 
gral. It is well known that if the moduli are quantised they are interpreted 
as differential operators in some representation of the universal enveloping 
algebra [35, 36, 76, 28, 29, 30, 31, 32]. Thus a solution of the loop equations 
can be found provided we find a trace (Q \ \ Q) that satisfies: 

0=({l\ ——q( x ,x;B) + 
+ / dy z B-\w-y)(Q | *(x,y;B) \Q)(Q\ *(y,x;B) | Q) (106) 

JC(x,x) 

for self-avoiding loops, that certainly happens if: 

= < fi I f$ < 107 > 

where we have introduced the " superscript to indicate that moduli are now 
differential operators on a Hilbert space. This is our new version of the mas- 
ter equation. In principle we could parameterise the state (Q |, that we are 
looking for, by means of the spectrum of the Hamiltonians Hi of the quan- 
tum Hitchin system (see appendix B), following the technique used in [28] 
to solve the Bethe ansatz equations for the quantum Hitchin Hamiltonians 
(on a punctured sphere). The basic idea there is to interpret Hitchin Hamil- 
tonians as elements of the centre of a Kac-Moody algebra at critical level, 
in a way that we recall briefly below (the trace property should be checked 
separately, using perhaps the quasi- conformal structure associated to the 
critical level [31, 32]). Langlands duality is a "far reaching generalisation of 
Artin-Hasse theory in number field" [77] that loosely speaking associates to 
every adelic automorphic representation on a local field a representation of 
its Galois group. Beilinson and Drinfeld started a program of understand- 
ing Langlands duality from a geometrical point of view. From this point 
of view to the Galois side corresponds a representation of the fundamental 
group of a Riemann surface, while to the automorphic side correspond cer- 
tain D-modules on the moduli stack of vector bundles on a Riemann surface. 
More precisely Beilinson and Drinfeld proved that the spectrum of the set of 
commuting operators obtained quantising the Hitchin Hamiltonians (i.e. the 
D-modules on the automorphic side) is parameterised by opers with values 
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in the dual Langlands group [26]. To the opers is then associated a repre- 
sentation of the fundamental group via the monodromy representation (this 
is the Galois side). The D-module associated to the quantisation of Hitchin 
bundles is in fact a representation of a Kac-Moody algebra of critical level, 
actually of its centre [28, 31, 32]. It is known that the common eigenstates 
of the set of Hitchin Hamiltonians can be obtained by the Bethe ansatz [28] 
(see also [76]). 

It is also clear that the first variation of the effective action computed at 
the master field can be considered as an element of the universal enveloping 
algebra of the global gauge group, because of translational invariance. But 
the universal enveloping algebra of any Lie group can be embedded into the 
enveloping algebra of the corresponding Kac-Moody algebra via the constant 
matrices. Therefore we can add to the Hitchin Hamiltonians, i.e. to the cen- 
tre of the universal enveloping algebra of the Kac-Moody algebra of critical 
level, just one operator in the universal enveloping algebra and still get a 
commutative family. We choose to add precisely the first variation of the 
effective action associated to the quantised Hitchin system. We can then 
solve the master equation: 



where (fl | varies over the family of the common eigenstates of the quantum 
Hitchin Hamiltonians that are parameterised by the moduli of opers for the 
dual Langlands group: 



Thus we come to the conclusion that while in sect. (4) the master field was 
parameterised by the moduli of the representations of the fundamental group 
with values in the global gauge group, here the vacuum is parameterised by 
the moduli of representations in the dual Langlands group. It is tempting to 
identify Langlands duality with 't Hooft electric-magnetic duality since the 
£7(1) subgroups of the gauge group are mapped into each other by Langlands 
duality precisely in the same way Abelian electric-magnetic duality does. We 
should mention that there is a natural physical interpretation of this construc- 
tion in terms of the conjectured existence of the gauge- field/string correspon- 
dence [33, 34]. The main difference, apart from the one already mentioned 



= {Q 




(108) 



Sfi(x) 



(Q | Ai = (ft | Hi 



(109) 
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about quantisation, is that if we insist to interpret Hitchin Hamiltonians as 
elements of the centre at critical level, there is no associated conformal ver- 
tex algebra, but rather only a quasi- conformal one into which the master 
field would live inside [31, 32]. We should also mention that the construction 
of sect. (4) furnishes, after quantisation, two different representations of the 
master field, related by a non-unitary similarity transformation. Indeed the 
master field b in the holomorphic gauge would be an element of a chiral vertex 
algebra, while the master field B in the unitary gauge would be an element 
of a non-chiral one. While this does not affect the Wilson loop, because of 
the required trace property of (Q \\ Q), some non-holomorphic information 
(i.e the harmonic metric) is necessary to compute the effective action T e ff. 



7 Conclusions. 



It is known that the von Neumann algebra generated by the master field of 
QCD in d dimensions, for d larger or equal to two, or by any equivalent 
Eguchi-Kawai reduction to any lower dimension, is a sub-algebra of a Cuntz 
algebra with more than one self-adjoint generator, that is not hyper-finite, 
i.e. it is not the weak limit of a sequence of finite dimensional matrix al- 
gebras. Yet functional techniques have shown that for QCD 2 some special 
observables, such as self-avoiding Wilson loops, admit a description in terms 
of a master field that generates a hyper-finite (Abelian) von Neumann alge- 
bra. We have shown in this paper that for planar self-avoiding Wilson loops, 
constructed by the connection B = (A z + nD u )dz + (A 2 — RD u )dz with k^O 
in a two dimensional Eguchi-Kawai model equivalent to the large- N limit of 
QCD4, the large- N loop equations admit as solution a master monodromy 
that, assuming existence of the large- N limit, generates a hyper-finite von 
Neumann algebra. An analogous construction holds for a partially twisted 
version, i.e. for the equivalent non-commutative gauge theory. In partic- 
ular the master algebra of the partially quenched or twisted EK model is 
the limit group algebra of a sequence of (possibly projective) representations 
(with values in the complexified gauge group) of the fundamental group of an 
arbitrarily punctured torus. As such it can be related to infinite dimensional 
Hitchin systems and in fact the loop equations written in terms of the " adelic 
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field of moduli" of the Hitchin systems: 



r 5T 
= / II D^' q DfI q ' exp(-r)(Tr(— V(x p , x p , B)) + 

+ / dy z d-\x p -y)Tr^{x p ,y-B))Tr^{y,x p -B))) (110) 

JC(x p ,Xp) 

are shown to be implied, for planar self-avoiding loops, by a master equation 
that determines the critical set of the following effective action: 

r e// = S YM - ^logDeti-A^B)) + logDet ^j^ + 



-t(C) l -jTr^d 2 x (111) 



Remarkably the first determinant in Eq.(lll) can be computed in a mani- 
festly four dimensional Euclidean (or Lorentz) invariant way and, once in- 
serted into the master equation, implies the correct value of the one-loop 
beta function up to terms of higher order in powers of the moduli fields. 
There is a dual construction, in the sense of the geometric Langlands cor- 
respondence, in which the Hitchin systems are quantised, the moduli are 
differential operators and the effective action is an element of the enveloping 
algebra of the global gauge group. Thus the loop equations: 

o = (n | — — m{x,x;B) | n) + 
oii(w) 

+ / dy z d- 1 (w-y)(Q\^(x,y;B)\Q)(Q\^(y,x;B)\Q) (112) 

JC(x,x) 

are transformed, for self-avoiding loops, into the condition that the master 
equation annihilates the vacuum: 

°-< n 'Sg (113) 

A possible physical interpretation of this quantised construction is in terms of 
a quasi-conformal (vertex) algebra relative to the conjectured string/gauge- 
field duality. 
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8 Appendix A 



In this appendix we follow [41]. 

A von Neumann algebra is a sub-algebra of the algebra of all bounded oper- 
ators on a Hilbert space that is closed under the adjoint operation (*) and 
under taking limits in the topology of weak convergence, the weak topol- 
ogy in short. The dual of a von Neumann algebra is the set of all bounded 
linear functionals over the algebra. It is usually assumed that the Hilbert 
space over which the algebra acts has a countable basis and that the von 
Neumann algebra has a separable dual. Von Neumann double commutant 
theorem asserts that a von Neumann algebra coincides with the commutant 
of its commutant. The commutant is the set of all bounded operators that 
commute with all the operators of the algebra. The centre of an algebra is 
the set of all the operators of the algebra that commute with all operators of 
the algebra. 

A von Neumann algebra is a factor if its centre is trivial, i.e. a multiple of the 
identity. Every von Neumann algebra decomposes as a direct integral of fac- 
tors. Two von Neumann algebras are algebraically isomorphic if there exists 
an invertible linear continuous map that preserves the "-algebraic structure. 
An element of the dual of a von Neumann algebra is usually called a state 
or a weight on the algebra. It is called positive if it maps positive operators 
(i.e operators whose spectrum lies in [0,oo[) into the positive numbers. A 
representation of a von Neumann algebra is a linear *-homomorphism into 
the algebra of bounded operators on a Hilbert space. Every positive weight, 
w, on a von Neumann algebra defines a representation of the algebra via the 
celebrated GNS representation. The GNS representation defines the scalar 
product on the Hilbert space of the representation as (A \ B) — w(A*B) and 
the action of the algebra on the Hilbert space in the obvious way. 
The basic result in the classification of von Neumann algebras is the exis- 
tence of an algebraic invariant of the algebra, the von Neumann dimension 
function, that associates to every orthogonal projector, P, of the algebra 
(P = P*,P = P 2 ) a non-negative number. Factors are classified according 
to the range of values that the dimension function may assume: 
type I: 0, 1, 2, ...oo 
type Ih: [0,1] 
type I loo: [0, oo] 
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type III: 0, oo . 

A von Neumann algebra is called finite if it is finite the range of its dimen- 
sion function. The theory of von Neumann algebras is a non-commutative 
far reaching generalisation of the integration theory and spectral theorem to 
which it reduces in the commutative case. Indeed every commutative von 
Neumann algebra is algebraically isomorphic to the algebra of all bounded 
measurable functions L°°(X,dfi) for some measure space X and thus every 
weight on a commutative von Neumann algebra is a measure. 
One of the most striking features of the von Neumann algebras is the occur- 
rence of a continuous dimension. In fact for type / and type II the dimension 
function can be identified with the trace (suitably normalised) of an orthog- 
onal projection. Therefore type / and II possess semi-finite traces, while 
type III has no non-trivial trace. A trace, r, is an element of the dual of 
the algebra that satisfies t(AB) = t(BA). If it exist it is unique for a factor 
up to a multiplicative constant. 

After years of efforts since their invention by Murray and von Neumann 
a general structure theorem has been achieved in the classification of the 
hyper-finite (or amenable or injective) factors with separable dual. A von 
Neumann algebra is hyper-finite if it is the closure of a sequence of finite 
dimensional matrix algebras. Though many von Neumann algebras satisfy 
this property, there are some interesting algebras that do not, as the Cuntz 
algebra with more than one self-adjoint generator. Here it is the complete 
list of the hyper-finite factors with separable dual. 
I n : the algebra of matrices of rank n. 

loo', the algebra of all bounded operators on a Hilbert space with countable 
basis. 

R: the unique hyper- finite factor of type II\. R has many equivalent real- 
izations: as the countable tensor product of n x n matrices equipped with 
the corresponding unique factorized normalised trace; as a tracial represen- 
tation of the canonical anti-commutation relations on an infinite dimensional 
Euclidean space; as a Krieger factor associated to an ergodic transformation 
with an invariant measure. 

i?o,i = R x loo'- the unique hyper-finite factor of type I loo- 
IIIq: Rw, the Krieger factor associated to an ergodic flow with a quasi- 
invariant measure. 

III\. R\ with A 7^ 0, 1, the Powers factors. 
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I Hi. Roo = R\i x R\ 2 with not rational, the algebra of canonical com- 
mutation relations in infinite dimensions. 

In particular the Powers construction of hyper-finite factors of all types in 
terms of pairs: (infinite inductive tensor product of matrix algebras, factor- 
ized state over the tensor product), whose type depends on certain conditions 
about the eigenvalue list of the states, turned out to be fundamental. 
We now discuss an application of the algebraic approach to the large- N limit 
that had some relevance for us in the introduction and in sect. (3). The 
simplest example of a hyper-finite master algebra is the one generated by 
functions of any single self-adjoint operator, that obviously is a commuta- 
tive one, equipped with an arbitrary tracial state. This may be the algebra 
of local observables of a one-matrix model or, less trivially, the sub-algebra 
generated by the master field for a self-avoiding loop of area A in the large 
TV-limit of QCD 2 : 

exp(iVA(a + a*)) (114) 

where a and a* are elements of the Cuntz algebra. In this commutative case 
hyper-finiteness follows from the Szego limit theorem for self- adjoint opera- 
tors, that depends essentially on the fact that Abelian C*-algebras are type 
I and thus every trace is uniformly approximately finite dimensional. In fact 
the interested reader may find in [18] a detailed study of how the approxi- 
mation properties of traces carry information about (tracial) representation 
theory of von Neumann algebras. 



9 Appendix B 

Hitchin has studied a non-canonical reduction of the four dimensional self- 
dual equations: 

F,u = (*F)^ (115) 

to two dimensions, requiring that the four dimensional gauge connection be 
constant along two of the four coordinates [5]. 

The reduction is non-canonical (i.e depending explicitly on a choice of the 
coordinate system) since Hitchin identifies the components of two polarisa- 
tions of the four-dimensional connection (A u , A a ) with the components of a 
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two dimensional one-form = + ip — A u dz + A a dz), the Higgs field, in 
order to write the reduced self-duality equations: 

F A - i^ 2 = 

8 A ifj = o 

d A i) = (116) 
as a flat condition for the non-Hermitian connection B = A + i^: 

F B = (117) 

plus a harmonic condition: 

d* A $> = (118) 

for the metric g associated to the one form ^ [78]: 

* = g- 1 dg (119) 

Hence the name of harmonic bundles for the Hitchin systems. The non- 
canonical nature of the Hitchin reduction matches nicely with the use that 
we have made of it in the partially quenched or twisted EK models in sect. (3). 
Hitchin has interpreted the reduced self-duality equations as the three mo- 
ment maps of a hyper-Kahler reduction for the Hamiltonian action of the 
local gauge group over the cotangent bundle of unitary connections (here we 
take the gauge group to be SU(N) or U(N)), labelled by the pairs (A,^). 
Hitchin equations can be generalised to include coadjoint orbits at points p 
of a divisor D: 

F A -im 2 = Y.^/ 2 \x-x P ) 

p 

d A ip = n P^ 2) ( x - x p) 
p 

d A ij = Yn p 6 (2) (x-x p ) (120) 
p 

It is known that for fixed eigenvalues of the Hermitian levels, the moduli 
space corresponding to taking the quotient by the action of the local compact 
gauge group is hyper-Kahler for n p nilpotent. Alternatively the moduli space 



47 



is hyper-Kahler for fixed coadjoint orbits, i.e. fixed monodromy around the 
punctures [69]. 

From our point of view, including coadjoint orbits on an arbitrary divisor is 
needed to provide a stratification of the functional integral that be, at the 
same time, dense in the sense of distributions and arising from geometry. In 
fact only for a collective field of the form: 

m° = E* (2) (^-^) 

P 

n = n p 8^ (x — x p ) 
p 

n = Y j n P 5 {2 \x-x p ) (121) 

p 

existence of solutions for the microcanonical ensemble of sect. (4) with a nice 
(i.e. Hausdorff) moduli space are guaranteed. 
The corresponding equation for B now reads: 

F B = /j,° + n-n (122) 

together with the harmonic condition: 

d* A ^J = n + n = v (123) 

The harmonic bundles defined by the moment map equations, with coad- 
joint orbits included, are called tame harmonic bundles since in this case the 
eigenvalues of ip are multi- valued holomorphic functions with poles at most 
of order one at the divisor of the coadjoint orbits. There are correspondences 
between tame harmonic bundles and several other objects of algebraic and 
topological origin [65, 66, 67]. Some of these correspondences play a funda- 
mental role in this paper because they allow us to give a precise meaning, on 
the dense set defined by Hitchin systems, to the formal change of variables 
that occur in the derivation of the master equation of sect. (5) from the loop 
equations. One way of understanding the existence of such correspondences 
comes from the following general result about hyper-Kahler quotients. As we 
have seen the hyper-Kahler quotient is obtained imposing the three moment 
maps and taking the quotient by the action of the compact gauge group. In 
this way we get the moduli space of tame harmonic bundles. However we can 
choose one of the three complex structures and the associated non-Hermitian 
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moment map and take the quotient under the action of the complexification 
of the local gauge group, getting the same moduli space. From Eq.(122) it 
is follows that to a tame harmonic bundle corresponds a flat connection and 
thus a local system, i.e. a representation of the fundamental group of the 
surface with a divisor. More precisely a local system is a fibre bundle with 
locally constant transition function and sections v given by solutions of the 
equations cIbv = for the flat connection B. The moduli of local systems 
can be written in terms of the monodromy of the connection B along the 
generators of the fundamental group of the Riemann surface with a divisor, 
subjected to the consistency condition: 

(a)tf (&)* (a) _1 tf (ft) -1 Y[V(D P ) = 1 (124) 

p 

after taking the quotient by the action of the global gauge group that acts by 
conjugation on the monodromies [75]. In addition to the flat connection it is 
possible to associate a flat holomorphic connection with regular singularities 
at the punctures: 

db' = iii (125) 

This is the holomorphic de Rham local system we refer to in the paper. 
Taking a reduction with respect to another complex structure we get the 
Dolbeault moduli space of Higgs bundles defined by: 

B A ip = n (126) 

The Dolbeault moduli do not play any role in sect. (5) but they will in sect. (6) 
when we consider a dual interpretation of the master equation based on 
the quantisation of Hitchin systems. In particular Hitchin Hamiltonians are 
polynomials in -0: 

^ = Tr(^) (127) 

For completeness we report a formula for the dimension of the moduli spaces 
of SU (N) Hitchin bundles on a Riemann surfaces with a divisor and nilpotent 
residue [79, 80, 81, 82, 83]. 

dimN = (2g - 2)(iV 2 - 1) + \ ]T f p (128) 

1 P eD 
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where f p = N 2 —J2 P fnt and m p the multiplicity of the weights of the parabolic 
flag at p [82]. We have seen in sect. (5) that for physical purposes, i.e. for 
calculating the string tension, it is more convenient to compute the unitary 
monodromy ty(x, x; A) on the master field rather than the B monodromy. We 
need therefore a way of separating the field A from \1/ in the loop equations. 
The natural way of doing it is to consider the new connection B K = A+i(nifj+ 
Rip) obtained by a suitable rescaling of the field tp. The A monodromy can 
then be recovered in principle in the limit n — > 0. We can thus compute the 
effective action for the microcanonical ensemble appropriate for B K and try to 
solve the master equation in the limit k — > 0. However we can do better. As 
it has been commented in the main text, it is in the spirit, if not in the letter 
of the master field idea, that the master field should be a k independent 
solution of the master equation up to gauge equivalence. Thus A and \1/ 
should be fixed points for the k action up to gauge equivalence. In finite 
dimensions it is known that such fixed points exist and are precisely those 
Hitchin bundles that arise as variations of Hodge structures [65, 66, 67, 68]. 
In addition for such finite dimensional bundles the monodromy of B K is k 
independent and coincides with ty(x,x; A) at the fixed point [68]. 



10 Appendix C 

In this appendix we study the loop equations for self-avoiding loops in large- 
iV QCD2 and we compare them with non-Abelian localisation. 
Localisation was discovered by Duistermaat and Heckmann [38] as a mean to 
compute exactly certain integrals by the saddle-point method despite they 
may not be Gaussian. 

It turned out that what was really involved in such integrals was some sort 
of cohomology. An infinite dimensional version of DH localisation due to 
Bismut [39, 40] was employed by Witten [37] for computing exactly the par- 
tition function and the cohomology ring of QCD 2 . 

For the sphere Witten localisation reduces to the saddle-point method. Since 
it depends on equivariant cohomology, Witten localisation applies to the par- 
tition function and to the observables belonging to the cohomology ring, but 
it does not apply to all the gauge invariant observables of QCD 2 . In par- 
ticular it does not apply to Wilson loops. Nevertheless it turns out that for 
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self-avoiding loops on a large sphere the result derived by large- N functional 
techniques and the one obtained were to hold localisation are surprisingly 
similar though they do not agree exactly. 

Indeed there is no reason why they should agree exactly since there is no 
reason why localisation should apply to Wilson loops. Here we use the tech- 
niques of sect. (3) and sect. (4) to explain where the discrepancy may arise. 
We have seen that our methods for loop equations lead to a result similar to 
localisation, in the sense that they transform loop equations for self-avoiding 
loops into a saddle-point problem for a certain effective action. Yet our 
saddle-point problem is in fact defined over a U(l) extension of the original 
theory, since as a consequence of the loop equations a central extension is 
present in the master equation. The functional techniques and two dimen- 
sional localisation instead would predict that there is no central extension at 
all. However we observed in the main text that, despite the master equa- 
tion is a priori defined over a central extension, we should require that the 
central extension vanishes due to a counter-term in the chiral Jacobian of 
the change of variables, as it actually happens provided we perform a finite 
conformal rescaling of the metric. In addition we showed that the master 
equation in our lattice or adelic version depends on the (regularized) value 
of the derivative §f (p) of the uniformization map at the base-point p of the 
loop, identified as a puncture of the bordered punctured surface internal to 
the loop. For QCD 2 our master equation coincides with the saddle-point 
equation obtained by localisation, up to the rescaling due to the mentioned 
derivative of the uniformization map. 

Thus it is precisely the existence of the uniformization map the source of the 
discrepancy between DH localisation in d = 2 and the present approach. In- 
deed we should remind the reader that the classical action is written in terms 
of fields living on the space-time manifold with global coordinates (z, z), while 
the functional determinants that arise in the effective action are defined in 
terms of fields living on the image of the uniformization map, with coordi- 
nates (t,t). Unfortunately we do not have a detailed information about the 
uniformization map, that by the way is singular at punctures, so that it must 
be regularized and suitably interpreted. What we do here is to find the value 
of §f (p) that is implied by the exact result obtained by functional methods. 
We find Iff (p)] 2 = Nc, where Nc is the number of punctures inside the loop. 
Localisation would imply instead |ff(p)| 2 = 1, somehow incorrectly unless 
the loop contains just one lattice point. 
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Let us present the detailed computation. For QCD2 it easy to follow the steps 
that lead to the master equation for QCD4. The microcanonical ensemble is 
described by the equation: 

F A -n = (129) 

that is the moment map at level \x for the Hamiltonian action of the local 
gauge group with respect to the symplectic structure: 



u) 



= J d 2 xTr(5A z A 5A Z ) (130) 



Alternatively Eq.(129) is the zero level condition for the gauge action on the 
gauge connection and the coadjoint orbits. Inserting the resolution of the 
identity: 

1 = J 8{F A -ii)Dii (131) 
into the functional integral: 

/N f 
DAD^i exp(-^ J TrF 2 A d 2 x)5(F A - (132) 

and changing variables in the functional integral from the moduli field /i of 
the unitary connections to the moduli // of the corresponding holomorphic 
local systems: 

dA' z -in' = (133) 
we get the effective action T for the moduli field //: 

= J £>//exp(-r) (134) 

where T is given by: 

r = — - / TrF\d 2 x - logDet'(d A A) + logDet'(^) (135) 

with the first determinant arising from localisation on a given level and the 
second determinant from the chiral change of variables (we assume that we 
are on a sphere or on a torus so that there are no non-local moduli). Here 
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the label ' for the functional determinants means that the zero modes coming 
from gauge invariance must be projected away. Then, as in sect. (4), the loop 
equations follow in the form: 

0= / Dfi'exp(-T)(Tr(—— -V(x,x;A)) + 
J d/i{w) 

+ [ dy x dr\w-y)Tr(V{x,y-,A))) (136) 

JC{x,x) 

From which we get the master equation: 

T e ff can be computed non-formally by adding the covariant gauge-fixing: 

d* A 5A = (138) 

with the corresponding Faddeev- Popov determinant and a choice of the met- 
ric: 

Teff = ^ J TrF 2 A d 2 x - hogDet(-A 1 (A)) + l -logDet{d\d A ) + 

+logDet(^) - e{C) 1 - f Trfid 2 x (139) 

l -S Lb Zi J 

However for computational purposes it is more convenient a different choice 
of gauge. If our Riemann surface is a sphere, the Hitchin system lives over 
a punctured sphere with punctures forming a uniform lattice (i.e. rotational 
invariant). The universal cover of a punctured sphere is the upper half plane 
H. Our punctured sphere is then conformally equivalent to a polygon in H 
with cusps on the rr-axis corresponding to the punctures. On H we choose 
the gauge A y = 0, that leaves, as a residual gauge symmetry, gauge transfor- 
mations that are y independent. In the gauge A y = the determinant due to 
localisation and the Faddeev-Popov determinant are both field independent 
and cancel each other. We can use the residual gauge symmetry to fix the 
gauge fi ch = (the label ch means the non-diagonal part) at the cusps (i.e 
the punctures) on the x-axis. Then the effective action reduces to: 

r e// = ^ E E K z)l2 + E E M/*f - h f 3 ) + 

+logDet{^ w ) - e (C)- J Trh^dH (140) 
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where the ^ and ^ superscript refer to the domain of definition of the lattice 
field, hp, and we have set a = with a the lattice spacing corresponding to 
the cutoff A of the theory, that comes from the product of delta functions at 
the same point in the classical action. The last two terms in Eq.(140) cancel 
each other (modulo 1) (up perhaps to terms of higher order in h®) provided 
the conformal factor in the metric is rescaled by a factor of it (see sect. (5)). 
The result coming from localisation would be instead: 

r = ^ E ^ 2 + E E io 9 (h^ - h^) (mi) 

that is similar up to the mentioned discrepancy in the field definitions. Should 
we compute the Wilson loop using the master field obtained by T rather than 
T e ff we would get the wrong result because of the following subtlety. The 
critical equation for T is just the saddle-point equation for the eigenvalues in 
the (z, z) coordinates. The master equation for T e ff is the critical equation 
for the lattice or adelic field of residues in the coordinates defined by the 
uniformization map, i.e. those labelled by t, after a conformal rescaling of 
a factor of tt to cancel the central term in the effective action. Let us sup- 
pose that we compute the Wilson loop according to the naive cohomological 
localisation: 

r N 

W(A) = Z- 1 / exp(--— E < - E log(fl ~ H)) E ex P(* E K) (142) 
Assuming translational invariance in the large- N limit we get: 

x^Texp(iA^) (143) 



W(A) = Z- 1 J exp(-^ E V 2 ~ E Nnlogih* - W)) x 



where N D is the number of punctures in our space-time manifold and N c 
the number of punctures internal to our loop C. By changing variables to 
9 % = N c h l we obtain: 

W(A) = Z-'J exp(-^^j £ 0- - g N D lo 9 (e' - IP)) x 

x^exp(^) (144) 
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The corresponding master equation for is: 

2A^=J^ < 145 > 

This equation should be compared with Eq.(16) that is the correct result 
obtained by functional methods: 

i^-^wha (146 » 

Since A = N c a 2 we see that Eq.(145) predicts a string tension wrong by a 
factor of Nc- 

In our master equation instead the functional derivative that acts on the 
classical term is rescaled by the factor of ||| (p)\ 2 while the functional deter- 
minants in the effective action are not. The master equation derived by T e ff 
now reads: 

m(p)\\ w = y ! (147) 



2\na 2 p " . Z i. h (t)i - h {t)j 



where we have included the conformal rescaling of a factor of n. We should 
however write it in terms of the the rotational invariant field h^: 

^IfMIV .. = y I ( 148 ) 

The equation for 9 l = Nc<h l p is then: 

"'gMIV =riL (14 9) 

that is the correct one provided |f|(p)| 2 = Nc- We can understand naively 
why |f;(p)| 2 should be of order of N c . Let us recall that our manifold has 
the topology of a punctured disk, with the boundary of the disk forming 
the " large horocycle" , while the punctures are regularized by means of small 
horocycles around them. Following [64, 84] , given a metric on a fatgraph with 
lengths k, there is a unique quadratic differential q such that the residue of 
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y/q at the puncture is k, the length of the horocycle around it. We will 
suppose that this picture extends to the case that the puncture is on the 
"large horocycle". We can identify the differential yfq with §f(p). After 
regularization it is of order of — that is proportional to the ratio between the 
lengths of the " large" and small horocycles. Its square is proportional to the 
ratio of the corresponding areas. But this is proportional to the number of 
punctures N c inside the "large horocycle". 



11 appendix D 

In this appendix we compute the divergent part of the localisation determi- 
nant. It is convenient to perform the computation in an indirect way, by 
means of a term by term comparison with the usual one-loop perturbative 
contribution to the effective action. For this purpose let us recall the one-loop 
perturbative result that is given by: 

r N N 

j DcIM /i exp(--Tr(c 2 )) e X p(--Tr(F%))6(D ll 5A ll - c)A FP = 

= Det^{-A A S^1 + 2iad F jDet(-A A ) (150) 

where we have inserted the gauge-fixing condition and the corresponding 
Faddeev- Popov determinant and by an abuse of language we have denoted 
with A the classical background field in the right hand side of Eq.(150). The 
perturbative computation of the one-loop beta function is the result of two 
contributions that are independent within logarithmic accuracy [57]. The 
orbital contribution gives rise to diamagnetism and to a positive term in the 
beta function: 

log(Det-*(-A A 8,»,l)Det(-A A )) = logDer\-A A ) = 

= T2(^ l0 ^I d ' xTr{F ^ < 151 » 

where it should be noticed the cancellation of two of the four polarisations 
between the first factor and the Faddev-Popov determinant. The spin con- 
tribution gives rise to paramagnetism and to a overwhelming negative term 
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in the beta function [57]: 



- A rr(i2ad F ^-^ A y l i2ad F ^X-/± A )- 1 ) 
12 N , / A 2 



log(-) f d 4 xTr{F, u f (152) 

Lb J 



12(2tt) 2 ^ v /i 

Hence the complete result is: 

11 N A 2 r A 2 r 

'Uj^y 109 ^ J dxTr ^ = / d 4 xTr{F, v f (153) 

In the case of the localisation determinant we have instead: 

i f 1 

DeH(-Ai) = limy DcDA^ exp(— -Tr(c 2 )) x 



1, 



x 



exp(--Tr((F" - //" ) 2 )<yp,M< - c)A FP (154) 

where we have inserted the gauge-fixing condition and the corresponding 
Faddeev-Popov determinant and by an abuse of language we have denoted 
by A the background field defined by the condition F~ v — n~ v = 0. The result 

for Det~z(— Ai) is thus: 

DeH(-Ai) = Der^{-A A 5^ U 1 + iad F -)Det(-A A ) (155) 

The coefficient of the spin contribution in this case is one half of the one in 
the perturbative case, but the spin term involves the anti-selfdual part of the 
field strength instead of the field strength itself. The orbital contribution is 
the same as in the perturbative case. Thus we have for the orbital part: 

1 AA A 2 r 1 
+ l2(2^ ,D9 ^>/ d V r < F ^ 156 > 

where we have used the decomposition: 



F»» = \f; u + \f; u (157) 
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For the spin part we get: 

^Tr(iad F - {-k A )- x iad F -{-k A )- x ) = 

= -T2W? l09{ $ ] I (158) 

We see that the spin part has the correct factor of 2 once it is written in 
terms of the anti-selfdual part of the field strength, a fact that we could have 
anticipated from start, to arrive at the final answer: 

that is the one reported is sect. (5). 
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